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Abstract. It has recently been shown that in every spatial dimension there exist 
precisely five distinct classes of topological insulators or superconductors. Within a 
given class, the different topological sectors can be distinguished, depending on the 
case, by a Z or a Z2 topological invariant. This is an exhaustive classification. Here 
we construct representatives of topological insulators and superconductors for all five 
classes and in arbitrary spatial dimension d, in terms of Dirac Hamiltonians. Using 
these representatives we demonstrate how topological insulators (superconductors) in 
different dimensions and different classes can be related via "dimensional reduction" by 
compactifying one or more spatial dimensions (in "Kaluza-Klein"-like fashion). For Z- 
topological insulators (superconductors) this proceeds by descending by one dimension 
at a time into a different class. The Z2-topological insulators (superconductors), 
on the other hand, are shown to be lower-dimensional descendants of parent Z- 
topological insulators in the same class, from which they inherit their topological 
properties. The 8-fold periodicity in dimension d that exists for topological insulators 
(superconductors) with Hamiltonians satisfying at least one reality condition (arising 
from time-reversal or charge-conjugation/particle-hole symmetries) is a reflection of 
the 8-fold periodicity of the spinor representations of the orthogonal groups SO(A'') 
(a form of Bott periodicity). Furthermore, we derive for general spatial dimensions 
a relation between the topological invariant that characterizes topological insulators 
and superconductors with chiral symmetry (i.e., the winding number) and the Chern- 
Simons invariant. For lower dimensional cases, this formula relates the winding number 
to the electric polarization (d = 1 spatial dimensions), or to the magnetoelectric 
polarizability (d = 3 spatial dimensions). Finally, we also discuss topological field 
theories describing the space time theory of linear responses in topological insulators 
(superconductors), and study how the presence of inversion symmetry modifies the 
classification of topological insulators (superconductors). 



PACS numbers: 73.43.-f, 74.20.Rp, 74.45. -fc, 72.25.Dc 
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1. Introduction 

Topological insulators (superconductors) are gapped phases of non-interacting fermions 
which exhibit topologically protected boundary modes. These boundary states are 
gapless, with extended wavefunctions, and protected against arbitrary deformations 
(or perturbations) of the Hamiltonian, as long as the generic symmetries (such as e.g. 
time-reversal symmetry) of the Hamiltonian are preserved and the bulk gap is not 
closed ||. The different phases within a given topological insulator (superconductor) 
are characterized by a topological invariant, which is either an integer Chern/winding 
number or a Z2 quantity. The integer quantum Hall effect (IQHE) is the best known 
example of such a phase. In the IQHE, the transverse (Hall) electrical conductance 
a^y carried by edge states is quantized and proportional to the topologically invariant 
Chern number As a consequence of the quantization of the Hall conductance axy, 

two different quantum Hall states with different values of a^y cannot be adiabatically 
connected without closing the energy gap in the bulk. Therefore they represent 
distinct "topologically ordered" phases that are separated by a quantum phase transition 
(figure |l]a) . 

Other examples of gapped topological phases are the (spinless) chiral px + Wy 
superconductor |^, which breaks time- reversal symmetry (T), and the quantum spin 
Hall (QSH) states [y, |10|, [Tl|, |T2|, |T3|, 0, |15[, which are time-reversal invariant. In the 



former case, the topological features in question are those of the fermionic quasiparticle 
excitations deep in the superconducting phase, whose dynamics is described by 
the Bogoliubov-de Gennes (BdG) Hamiltonian I]. The BdG Hamiltonian of any 
superconductor has, by construction, a "built-in" charge-conjugation (or: particle- 
hole) symmetry (C). In analogy to the IQHE, the different topological phases of the 



Px + iPy superconductor can be labeled by an integer (Z) topological invariant [16 
The corresponding quantized conductance is the transverse (Hall) thermal, or "Leduc- 
Righi" conductance (divided by temperature). The quantum spin Hall effect (QSHE), 
on the other hand, which occurs in two-dimensional and three-dimensional time-reversal 



invariant insulators, is characterized by a Z2 topological number [y, [T^, |T3|, |TJ] . This 
binary quantity distinguishes the non-trivial state (z/q = 1), whose boundary modes 
consist of an odd number of Dirac fermion modes (i.e., odd number of Kramers' pairs), 
from the trivial state (z/q = 0), which is characterized by an even number of Kramers' 
pairs of boundary (surface or edge) states (figure |l]b). A physical realization of the QSHE 



was theoretically predicted and experimentally observed [|T8], in HgTe/(Hg,Ce)Te 

X In the presence of disorder (= lack of translational symmetry) a sufficient condition is that all bulk 
states near the Fermi energy are localized Q. 

§ A profound implication of a non-vanishing Chern number is that exponentially localized Wannier 
wave functions cannot be constructed [§[ |[ D- Yet another manifestation of the non-trivial nature 
of the band structure when a^y 7^ appears in the scaling of the entanglement entropy (and the 
entanglement entropy spectrum) 

II Physically, the quasiparticles are responsible for heat (energy) transport, and not for the electrical 
transport properties of the superconductor. 
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(b) 



quantum phase transition 



quantum phase transition 



Figure 1. Topological distinction among quantum ground states, (a) Z classification 
and (b) Z2 classification. 



semiconductor quantum wells. Subsequently, based both on theoretical considerations 



||15| , pO|] and experimental measurements ^ the three-dimensional Z2 

topological insulator phase was shown to occur in Bismuth-related materials, such as 
BiSb alloys, Bi2Se3 and Bi2Te3 . 

The key difference among the IQHE, the QSHE and three-dimensional Z2 
topological insulator, and the chiral Px + iPj/ superconductor lies in their generic 
symmetry properties (to be discussed in more detail below). All these states are 
topological in the sense that they cannot be continuously deformed into a trivial 
insulating state (i.e., a state without any gapless boundary modes) while keeping 
the generic symmetries of the system intact and without closing the bulk energy 
gap. The key property of the QSH states, for example, is the spin-orbit interaction 
in combination with time-reversal symmetry (T). The latter protects the boundary 
states from acquiring a gap or becoming localized in the presence of disorder (= terms 
added to the Hamiltonian which break translational symmetry). The chiral px + ipy 
superconductor, on the other hand, is characterized by C but the lack of T. 

The aforementioned examples of "symmetry protected" topological states are 
part of a larger scheme|5| that fully classifies ^ all topological insulators 

(superconductors) in terms of symmetry and spatial dimension. This classification 
scheme is summarized in table ^. The first column in this table provides the complete 
list of all possible "symmetry classes" of single-particle Hamiltonians. There are 
precisely ten such "symmetry classes", a fundamental result due to Zirnbauer, and 



Altland and Zirnbauer ||29|, These authors recognized that there is a one- 

to-one correspondence between single-particle Hamiltonians and the set of ("large") 
symmetric spaces, of which there are precisely ten, a classic result obtained in 1926 by 

^ A complete classification of all topological insulators (superconductors) in dimensionalities up 
to d = 3 was given in Ref. A systematic regularity (periodicity) of the classification as the 



dimensionality is varied was discovered by Kitaev in Ref. [27 for all dimensions through the use of 
K-Theory. A certain systematic pattern as the dimensionality is varied was discovered for some of the 
topological insulators by Qi, Hughes and Zhang in Ref. using a Chern-Simons action in extended 
space. Below we will use a version of the ideas of "dimensional reduction" employed in Ref. pa|. 
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the mathematician Ehe Cartan. The work of Altland and Zirnbauer extends, indeed 
completes, the much earher work by Wigner and Dyson most well known from the 
context of random matrix theory as the three "unitary, orthogonal, and symplectic" 
symmetry classes |]32|. Before discussing in more detail the result of the classification 
in table 0, we will first briefly explain the notion of "symmetry class" and the ten-fold 
list of Hamiltonians, which is the framework within which the classification scheme in 
table ^ is formulated. The reader familiar with this topic may skip this subsection. 



1.1. Review: Classification of generic Hamiltonians - "The Ten-Fold Way" 

Consider the gapped first-quantized Hamiltonian describing the topological insulator 
(superconductor) in the bulk {d spatial dimensions). The topological features of 
this Hamiltonian, which we are interested in characterizing, are (by definition) not 
changed if we modify (deform) the Hamiltonian by adding terms (or perturbations) to 
it which break any translational symmetry that may be present. Thus the topological 
features must be properties of a translationally not invariant Hamiltonian. Seeking 
a classification scheme, one must have some framework within which to classify such 
translationally not invariant Hamiltonians. Such a classification cannot involve the 
notion of ordinary symmetries, i.e., of unitary operators which commute with the 
Hamiltonian. Consider, e.g., (Pauli-) spin rotation symmetry, ubiquitous in condensed 
matter systems. The notion of spin-rotation invariance can be eliminated by writing the 
Hamiltonian in block form, and by focussing attention on a block of the Hamiltonian. 
Such a block decomposition can be performed for any symmetry which commutes with 
the Hamiltonian. The notion of "symmetry classes" mentioned above refers to the 
properties of these blocks of symmetry- less "irreducible" Hamiltonians: as it turns out, 
there are only ten of them. The basic idea why this is the case is simple to understand. 
The only properties that the blocks can satisfy are certain reality conditions that follow 
from the "extremely generic symmetries" of time-reversal and charge- conjugation (or: 
particle-hole symmetry). These are not symmetries in the above sense, since they are 
both represented by aniz-unitary operators, when acting on the single-particle Hilbert 
space. Specifically, consider a general system of non-interacting fermions described by a 
"second-quantized" Hamiltonian H . For a non-superconducting system, e.g., this reads 

^ = E^A^AB^B, (1) 

A,B 

with fermion creation and annihilation operators satisfying canonical anti-commutation 
relations 

{^^,4} = 5^,5. (2) 

Here we imagine for convenience of notation that we have "regularized" the system 
on a lattice, and A.,B are combined labels for the lattice sites and if relevant, of 
additional quantum numbers such as e.g., a Pauli-spin quantum number (e.g., A = (i, a) 
with a,b = ±1/2). Then 'Ha,b is a.n N x N matrix, the "first quantized" Hamiltonian. 
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(Similarly, a superconducting system is described by a BdG Hamiltonian for which we 
use the Nambu-spinor instead of complex fermion operators ^Z'^, and whose first 
quantized form is again a matrix "H when discretized on a lattice.) 

Now, time reversal symmetry can be expressed in terms of "H: the system is invariant 
under time reversal symmetry if and only if the complex conjugate of the first quantized 
Hamiltonian l-L* is equal to 7/ up to a unitary rotation Ut, i.e. 

r : U^H* Ut = +n. (3) 

Moreover, the system is invariant under charge- conjugation (or: particle-hole) symmetry 
if and only if the complex conjugate of the Hamiltonian "H* = is equal to minus l-L 
up to a unitary rotation f/c, i.e. 

C: UlU* Uc = -H. (4) 



(This property may be less familiar, but it is easy to check that it is a 
characterization of charge-conjugation (particle-hole) symmetry for non-interacting 
systems of fermions). A look at equations (0,§) reveals that T and C, when acting 
on the single particle Hilbert space, are not unitary symmetries, but rather reality 
conditions on the Hamiltonian "H modulo unitary rotations Q Q. 

In second-quantized language, time-reversal and particle-hole operations can be written in terms of 
their action on the canonical fermion creation and annihilation operators, 

T^aT-i = iUT)A.B ^B, C^aC-^ = (Uc)a,b ^l- (5) 

B B 

While the particle-hole transformation is unitary, the time-reversal operation is anti- unitary, TW^^ = 
— i. The system is time- reversal invariant (particle- hole symmetric), if and only if THT^^ = H 
{CHC~^ = H). This leads directly to the conditions (||) and (^) for the first quantized Hamiltonian. 
Note that THT-^ = H implies T^A{t)T-^ = Te+'^'V^e-'^'r-i = Es {Ut)a,b 4'Bi-t). 
Iterating T and C twice, one obtains T^i/jaT^^ = J2b {UtUt)a,b "05, and C^VaC^^ = 
TliB {U^c^c)a,b "ipB- When acting on the first quantized Hamiltonian this reads {U^Ut)'^T-L{U^Ut)= 'H 
and {UqUc)'^T-L{UqUc)= H, respectively. The first quantized Hamiltonians H are seen (below) to run 
over an irreducible representation space, and thus (U^Ut) and {U^,Uc) are both multiples of the 
identity matrix In (by Schur's lemma). Since Ut and Uc arc unitary matrices, there are only two 
possibilities for each, i.e., U^Ut = ^In and UqUc = ^In- The time-reversal operation T and the 
particlc-holc transformation C can then each square to plus or to minus the identity, — ±1, and 
C2 =±1. 

* It may also be worth noting that we may assume without loss of generality that there is only a 
single time-reversal operator T and a single charge-conjugation operator C. If the (first quantized) 
Hamiltonian H was invariant under, say, two charge-conjugation operations Ci and C2, then the 
composition Ci ■ C2 of these two symmetry operations would be a unitary symmetry when acting on 
the first quantized Hamiltonian "H, i.e., the product Uci ■ Uq^ would commute with Ji. By bringing the 
Hamiltonian H in block form, Uci ■ U^^ would then be constant on each block. Thus, on each block Uc\ 
and would then be trivially related to each other, and it would suffice to consider one of the two 
charge conjugation operations. - On the other hand, note that the product T-C corresponds to a unitary 
symmetry operation when acting on the first quantized Hamiltonian H. But in this case the unitary 
matrix Ut ■ does not commute but anti-commutcs with H. Therefore, T ■ C does not correspond to 
an an "ordinary" symmetry otH. It is for this reason that we need to consider the product T-C (called 
"chiral" , or "sublattice" symmetry S below) as an additional essential ingredient for the classification 
of the blocks, besides time- reversal T and charge-conjugation (particle- hole) symmetries C. 
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Now it is easy to see that there are only ten possible ways for a system to respond to 
time-reversal and charge-conjugation (particle-hole symmetry) operations. As for time- 
reversal symmetry (T), the Hamiltonian can either be (i): not time-reversal invariant, 
in which case we write T = 0, or (ii): it may be time- reversal invariant and the anti- 
unitary time-reversal symmetry operator T squares to plus the identity operator, in 
which case we write T = +1, or (iii): it may be time-reversal invariant and the anti- 
unitary time-reversal symmetry operator T squares to minus the identity, in which case 
we write T = — 1. Similarly, there are three possible ways for the Hamiltonian "H to 
respond to charge- conjugation (particle-hole symmetry) C (again, C may square to plus 
or minus the identity operator). For these three possibilities we write C = 0, +1,— 1. 
Hence, there are 3x3 = 9 possible ways for a Hamiltonian % to respond to time- 
reversal and charge-conjugation (particle-hole transformation). These are not yet all 
ten cases because it is also necessary to consider the behavior of the Hamiltonian under 
the product S = T ■ C which is a unitary operation. A moment's thought (see table |l|) 
shows that for 8 of the 9 possibilities the behavior of the Hamiltonian under the product 
S = T ■ C is uniquely fixed ||. (We write S = if the operation S is not a symmetry 
of the Hamiltonian, and S = 1 if it is.) The only case when the behavior under the 
combined transformation S = T ■ C is not determined by the behavior under T and C 
is the case where T = and C = 0. In this case either S = or S = 1 is possible. This 
then yields (3x3 — 1) + 2 = 10 possible behaviors of the Hamiltonian. 

The list of ten possible behaviors of the first quantized Hamiltonian under T, C and 
S is listed in table ^ These are the ten generic symmetry classes (the "ten-fold way" ) 
which are the framework within which the classification scheme of topological insulators 
(superconductors) is formulated. 

Let us first point out a very general structure seen in table |T]. This is listed 
in the column entitled "Hamiltonian". When the first quantized Hamiltonian "H is 
"regularized" (or: "put on") a finite lattice, it becomes an N x N matrix (as discussed 
above). The entries in the column "Hamiltonian" specify the type of N x N matrix that 
the quantum mechanical time-evolution operator exp{itl-L) is. For example, for systems 
which have no time- reversal or charge- conjugation symmetry properties at all, i.e., for 
which T = 0, C = 0, S = 0, which are listed in the first row of the table, there are no 
constraints on the Hamiltonian except for Hermiticity. Thus "H is a generic Hermitian 
matrix and the time-evolution operator is a generic unitary matrix, so that exp(it'H) is 
an element of the unitary group U(A^) of unitary N x N matrices. By imposing time- 
reversal symmetry (for a system that has, e.g., no other degree of freedom such as, e.g., 
spin), there exists a basis in which "H is represented by a real symmetric N x N matrix. 
This, in turn, can be expressed as saying that the time-evolution operator is an element 
of the coset of groups, U(A^)/0(A^). All other entries of the column "Hamiltonian" can 

ft The symmetry operation S is sometimes called "sublatticc symmetry" , hence the notation S. 
However, in many instances S is not realized as a sublattice symmetry, but is just simply the product 
of T and C. Therefore the term "chiral symmetry" to describe the symmetry operation S is sometimes 
more appropriate. 
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Cartan label 


T 


c 


s 


Hamiltonian 


G/H (ferm. NLctM) 


A (unitary) 











U(iV) 


U(2n)/U(n) x U(n) 


AI (orthogonal) 


+1 








V[N)/0{N) 


Sp(2n)/Sp(n) x Sp(7i) 


All (symplectic) 


-1 








U(2Af)/Sp(2Af) 


0(2n)/0(n) x 0(n) 


AIII (ch. unit.) 








1 


\]{N + M)/\]{N) X U(Af) 


U(n) 


BDI (ch. orth.) 


+1 


+1 


1 


0{N + M)/0{N) X 0(M) 


U(2n)/Sp(2n) 


CII (ch. sympl.) 


-1 


-1 


1 


Sp(iV + A/)/Sp(iV) X Sp(M) 


U(2n)/0(2n) 


D (BdG) 





+1 





S0(2iV) 


0(2n)/U(n) 


C (BdG) 





-1 





Sp(2iV) 


Sp(2n)/U(n) 


Dili (BdG) 


-1 


+1 


1 


S0(2iV)/U(A^) 


0(2n) 


CI (BdG) 


+1 


-1 


1 


Sp(2iV)/U(A^) 


Sp(2n) 



Table 1. Listed are the ten generic symmetry classes of single-particle Hamiltonians 
classified according to their behavior under time- reversal symmetry (7~), charge- 
conjugation (or: particle-hole) symmetry (C), as well as "sublattice" (or: "chiral") 
symmetry [S). The labels T, C and S, represent the presence/absence of time- 
reversal, particle-hole, and chiral symmetries, respectively, as well as the types of these 
symmetries. The column entitled "Hamiltonian" lists, for each of the ten symmetry 
classes, the symmetric space of which the quantum mechanical time-evolution operator 
exp(it7^) is an element. The column "Cartan label" is the name given to the 
corresponding symmetric space listed in the column "Hamiltonian" in Elie Cartan's 
classification scheme (dating back to the year 1926). The last column entitled "G/7J 
(ferm. NLcrM)" lists the (compact sectors of the) target space of the NLuM describing 
Anderson localization physics at long wavelength in this given symmetry class. 

be obtained from analogous considerations ||||. What is interesting about this column 
is that its entries run precisely over what is known as the complete set of ten ( "large" ) 
symmetric spaces U|, classified in 1926 in fundamental work by the mathematician Elie 
Cartan. Thus, as the first quantized Hamiltonian runs over all ten possible symmetry 
classes, the corresponding quantum mechanical time-evolution operator runs over all ten 
symmetric spaces. Thus, the appearance of the Cartan symmetric spaces is a reflection 
of fundamental aspects of (single-particle) quantum mechanics. We will discuss the last 
column entitled "G/if (ferm. NLaM)" in the following subsection. 

ft Possible realizations of the chiral symmetry classes AIII, BDI, CII possessing time-evolution 
operators in table ^ with N ^ M are tight-binding models on bipartite graphs whose two (disjoint) 
subgraphs contain and M lattice sites. 

f A symmetric space is a finite-dimensional Riemannian manifold of constant curvature (its Riemann 
curvature tensor is covariantly constant) which has only one parameter, its radius of curvature. There 
are also so-called exceptional symmetric spaces which, however, are not relevant for the problem at 
hand, because for them the number N would be a fixed finite number, which would prevent us from 
being able to take the thermodynamic (infinite- volume) limit of interest for all the physical systems 
under consideration. 



CONTENTS 



10 



1.2. Review: Classification of Topological Insulators (Superconductors) 



The approach used in Ref. to classify all possible topological insulators 

(superconductors) rested on the existence of protected extended degrees of freedom 
at the system's boundary. These are, in particular, also protected in the presence of 
arbitrarily strong perturbations at the boundary which break translational symmetry 
(commonly referred to as "random", or "disordered"). The existence of extended, 
gapless degrees of freedom even in strongly random fermionic systems is highly unusual, 
because of the phenomenon of Anderson localization |||. Thus, the degrees of freedom at 
the boundary of topological insulators (superconductors) must be of a very special kind, 
in that they entirely evade the phenomenon of Anderson localization. Our approach 
consists in classifying precisely all those problems of (non-interacting) fermionic systems 
which completely evade the phenomenon of Anderson localization. We have completed 
this task in Ref. ||2^. Thus, we have reduced the problem of classifying all topological 
insulators (superconductors) in d spatial bulk dimensions to a classification problem of 
Anderson localization at the [d — 1) dimensional boundary. By solving the mentioned 
problem of Anderson localization, we have thereby solved the problem of classifying all 
topological insulators (superconductors). 

We will now focus attention on the above-mentioned problem of Anderson 
localization at the {d — 1) dimensional boundary of a d dimensional topological insulator 
(superconductor). Specifically, we will now review a solution of this problem pB|, 



which allows us to see directly the dependence on dimensionality d of the classification of 
topological insulators (superconductors). In the later, main part of this article, we will 
present another point of view of this dependence on dimensionality d (using "dimensional 
reduction"). 

The theoretical description of problems of Anderson localization is well known to 
be very systematic and geometrical A problem of Anderson localization is in 

general described by a random Hamiltonian (i.e., one that lacks translational symmetry). 
That Hamiltonian will be in one of the ten symmetry classes listed in table |l| and 
we are currently focusing on Hamiltonians describing the boundary of the topological 
insulator (superconductor). Now, as it turns out, at long length scales (much larger than 
the "mean free path") a description in terms of a "non-linear-sigma-model" (NLcrM) 
emerges. A NLcrM is a system like that describing the classical statistical mechanics 
of a Heisenberg magnet. The only difference is that while the magnet is formulated in 
terms of unit vector spins, pointing to the surface of a two-dimensional sphere, for a 
general NLaM that spin is replaced || by an element of one of the ten symmetric spaces 
listed in the last column of table called the "target space" of the NLcrM, denoted by 

X This is the phenomenon that, at least for sufficiently strong disorder potentials, spatially extended 
{— delocalized) eigenstates of the Hamiltonian tend to become localized (i.e., exponentially decaying 
in space) ||, ||, ||, ||. 

§ The two-dimensional sphere S"^ is a particularly simple example of a symmetric space, namely it can 
be written as the space S"^ = U(2)/U(l) x U(l), i.e., the first row in the last column of table 0, when 
n = 1. 
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G/H 1]. For a given symmetry class of the original Hamiltonian, whose time-evolution 
operator is characterized by the penultimate column, the specific "target space" that 
needs to be used for the NLaM describing Anderson localization in this symmetry class 
is listed in the last column of table ^ (see also appendix 

Now, the NLcrM on the {d — 1) dimensional boundary of the d dimensional 
topological insulator (superconductor) completely evades Anderson localization if a 
certain extra term of topological origin can be added to the action of the NLcrM which 
has no adjustable parameter ||26|, Whether such an extra term is allowed depends 
on (i): the "target space" of the NLcrM in the symmetry class in question^ and (ii): 
the dimensionality d := {d — 1) oi the boundary on which the NLcrM is defined. There 
are only three terms of topological origin which can possibly be added to the action of 
the NLcrM: these are a ^-term (Pruisken term) a Z2 topological term 



or a 



Wess-Zumino-Witten (WZW) term It is the homotopy groups of the NLcrM 

target spaces G/H which determine whether it is possible to add such a topological term 
to a given NLcrM action (see table H). Specifically, a 6'-term (Pruisken term) can appear 
when T[g{G/H) = Z. Similarly, a Z2 topological term is allowed when n^lG/H) = Z2, 
and a WZW term can be included Q when 7rj_^_i{G/H) = Z. Note that, on the one 
hand, one obtains, upon addition of a ^-term (Pruisken term), a one-parameter family 
of theories (on the boundary) depending on the value of 6, all of which reside in the 
same symmetry class. On the other hand, however, it is known that only for a special 
value of the parameter 6 Anderson localization is avoided; for generic values of 6, this 
is not the case. It is for this reason that the ability to add a ^-term (Pruisken term) 
is not of interest to the question we are asking. One is therefore left only with a Z2 
topological term or a WZW term as the only terms of topological origin which have no 
freely adjustable parameter, and which are thus of relevance here. 

The homotopy groups for all ten NLcrM target spaces G/H listed in the last column 
of table |l| are well known from the literature (see, e.g., Ref. [^] for a summary) and 
this information is summarized in table for the convenience of the reader. In order 
to make a certain regular structure of this table apparent, the rows of table have 
been re-ordered in a specific way Part of this re-ordering is a subdivision of all 



symmetry classes into what is called ^^complex case" and "rea/ case" in table |^. The 
physical origin of this subdivision is simple to understand. In the category ^^complex 
case" appear precisely those symmetry classes in which there is no reality condition 
(T or C) whatsoever imposed on the Hamiltonian. We see from table |l| that these are 
the classes which carry Cartan label A and AIII. The Hamiltonians in these symmetry 
classes are therefore "complex" (in this sense). The other category called "rea/ case" in 

II Both G and H are classical Lie groups, and iJ is a maximal subgroup of G. 

% The relevant "target space" is, as already mentioned above, listed in the last column of table |l|. 

Observe that, while the homotopy group determines if a term of topological origin for a given NLctM 
action is allowed in principle, it depends on the specific disorder model considered, whether such a term 
is actually present. In any case, if a term of topological origin is possible in a specific dimension and 
symmetry class, then a corresponding topological insulator (superconductor) can exist in this symmetry 
class (in one dimension higher). 
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table 1^ consists precisely of all the other eight symmetry classes which have the property 
that there is at least one reality condition [T or C - cf. equations (0^)] imposed on the 
Hamiltonian (see table In this sense the Hamiltonians in these symmetry classes are 
therefore "real". 

complex case: 

AZ G/H J= d = l d = 2 d = 3 d = 4 d = 5 d = 6 d = 7 

A V{N+M)/ViN)xV{M) 4- Z < — Z < — Z < — Z ^ 

AIII \j{N) ^ — Z ^ — Z ^ — Z ^ — Z 



real case: 



AZ G/H d = d=l d = 2 d = 3 d = 4 d = 5 d = 6 d = 7 

Al Sp{N+M)/Sp(N)xSp(M) ^ Z < — Z Z2 Z2 ^ 

BDI u(2iv)/sp(2iv) ^ — Z ^ — Z Z^ Z^ 

~15 o(27V)/u(iv) Z^ ^ — Z ^ — Z Z^ 

Dili o{iv^ % % ^ — Z ^ — Z 

All 0{N+M)/0(N)xO{M) ^ Z Z2 Z2 i — Z ^ 

CII u(7V)/o(JV) ^ — Z Z^ Z^ ^ — Z 

~~C Sp(27V)/u(7V) ^ — Z Z^ Z^ ^ — Z 

~CI Sp(27V) ^ — Z Z^ Z^ ^ — Z 



Table 2. Tabic of homotopy groups tt^{G/H) for symmetric spaces G/H, taken from 
the standard mathematical literature (see e.g. Ref. |Q for a summary). [Here, N must 
be sufficiently large for a given d. Cartan labels of those symmetry classes invariant 
under the chiral symmetry operation S = T ■ C from table |l| are indicated by bold 
face letters.]. The pattern continues for higher d, with periodicity 2 for the complex 
case, and 8 for the real case. The entries corresponding to Z topological insulators 
(superconductors) in (d + 1) dimensions (see e.g. table ^) are indicated by blue color 
boldface symbols. There is a (d + l)-dimensional Z2 topological insulator, whenever 
TT^{G/H) = Z2 (also indicated in blue). A {d + l)-dimensional Z topological insulator, 
on the other hand, can be realized, whenever 7rj_|_]^ (G/iJ) = Z. - This is a way |2^, ^ 
to directly relate the classification of topological insulators (superconductors) to the 
table of homotopy groups. 

Table ^ now tells us directly in which symmetry class there exist topological 
insulators or superconductors (and of which type, Z2 or Z): according to the above 
discussion we know that a topological insulator (superconductor) exists in a given 
symmetry class in d = d + 1 spatial dimensions if and only if the target space of the 
NLcxM on the d-dimensional boundary allows either (a): for Z2 topological term, which 
is the case when tt^{G/H) = 7Td-i{G/H) = Z2, or (b): allows for a WZW term, which is 
the case when ^^{G/ H) = -K^^iiG/ H) = Z. By using this rule in conjunction with table 
of homotopy groups, we arrive with the help of table ^ at the table ^ of topological 
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Table 3. Classification of topological insulators and superconductors as a function of 
spatial dimension d and symmetry class, indicated by the "Cartan label" (first column). 
The definition of the ten generic symmetry classes of single particle Hamiltonians 
(due to Altland and Zirnbauer[^ ^) is given in table ^ The symmetry classes 
are grouped in two separate lists, complex and real cases, depending on whether the 
Hamiltonian is complex, or whether one (or more) reality conditions (arising from time- 
reversal or charge-conjugation symmetries) are imposed on it; the symmetry classes 
are ordered in such a way that a periodic pattern in dimensionality becomes visible 
[ p7| . (See also the discussion in subsection 1.1 and table ^.) The symbols Z and 
Z2 indicate that the topologically distinct phases within a given symmetry class of 
topological insulators (superconductors) are characterized by an integer invariant (Z), 
or a Z2 quantity, respectively. The symbol "0" denotes the case when there exists 
no topological insulator (superconductor), i.e., when all quantum ground states are 
topologically equivalent to the trivial state. 



insulators and superconductors Q. 

A look at table |^ reveals that in each spatial dimension there exist five distinct 
classes of topological insulators (superconductors), three of which are characterized by 
an integral (Z) topological number, while the remaining two possess a binary (Z2) 
topological quantity |[ 

The topological insulators (superconductors) appear in table ^ along diagonal lines. 
That is, as the spatial dimension is increased by one, locations where the topological 

* To be explicit, one has to move all the entries Z in tabic || into the locations indicated by the arrows, 
and has to replace the column label dhy d = d-\-l. The result is table |[ 

tt Note that while d = 0, 1, 2, 3-dimensional systems are of direct physical relevance, higher dimensional 
topological states might be of interest indirectly, because, for example, some of the additional 
components of momentum in a higher dimensional space may be interpreted as adiabatic parameters 
(external parameters on which the Hamiltonian depends, and which can be changed adiabatically, 
traversing closed paths in parameter space - sometimes referred to as adiabatic "pumping processes" ) . 
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insulators (superconductors) appear in the table shift down by one column. This entire 
pattern descends directly from an analogous one appearing in the table ^ of homotopy 
groups. 

In the main part of the present paper we will give another derivation of this 
periodicity structure of table 0, by "dimensional reduction" and by considering the 
Dirac Hamiltonian representatives in each of the symmetry classes in which there exist 
topological insulators (superconductors) in a given spatial dimension. By running the 
argument backwards, this can then, in turn, also be viewed as another route to obtain 
the table || of homotopy groups. 

Let us pause for a moment to point out a few interesting aspects relating to table 
First, it is interesting to note, that the periodic structure of table ^ is quite analogous 
to the dimensional hierarchy of anomalies ( "anomaly ladder" ) known in gauge theories, 
e.g., the (3+l)-dimensional Abelian anomaly, the (2+l)-dimensional parity anomaly, 
and the (H-l)-dimensional non- Abelian anomaly This may not appear to 

be completely unexpected, considering that the IQHE is known as a condensed matter 
realization of the parity "anomaly" |f|[ 

Second, we want to draw the reader's attention to the fact that for a given symmetry 
class in table |^ any Z2 topological insulator (superconductor) always appears as part of 
a "triplet" that consists of a rf-dimensional Z topological insulator (superconductor) and 
two Z2 topological insulators (superconductors) in {d — 1) and {d — 2) dimensions. This 
suggests that the topological characteristics of the members of such a "triplet" are closely 
related and indeed, it was shown in Ref. that the Z2 classification in symmetry class 
All m d = 2 and 3 can be derived from the four-dimensional Z topological insulator 
by a process of "dimensional reduction" . The same procedure has been applied ||28| to 
derive Z2 classifications in those symmetry classes that do not possess a form of chiral 
symmetry (i.e., for which S = 0) - there are five of them. In section ^ we will extend 
this approach to the Z2 topological insulators (superconductors) in the remaining five 
symmetry classes which possess a form of chiral symmetry (i.e., for which S = 1). 

Let us now discuss explicit forms of the "terms of topological origin" which can 
be added to the action of the NLaM with "target space" G/H , at the boundary of the 
topological insulator (superconductor). 



1.2.1. Terms of WZW type. Here there is a pattern, alternating in the dimension 
d of the boundary, for the way these topological terms can be constructed, (i): For 
a. d = odd-dimensional boundary, and for the NLaMs which describe the Anderson- 
delocalized boundary of the topological bulk, the field (denoted by Q G G/H below |f] 

ff The IQHE is not quite an anomaly in the sense that the parity and T is exphcitly broken. 
However, the quantization of Uxy in the IQHE can be viewed as essentially the same phenomenon 
as the appearance of the Chern-Simons term in QED in Z? = 2 + 1-dimensional space-time, where 
massless Dirac fermions are coupled with the electromagnetic U(l) gauge field. 

f While looking similar, the field Q in the NLcrM has nothing to do with the spectral projector in 
momentum space defined in equation (M). 
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- cf. last column of table of the NLaM field theory is a Hermitean matrix field 
which satisfies certain constraints. The WZW term responsible for the lack of Anderson 
localization takes on the form of the (i- dimensional integral over the boundary of the 
(i- dimensional Chern-Simons form (compare equation ( p2D of the main text below). 
Alternatively, it can be written also as an integral over a ci + 1-dimensional region whose 
boundary coincides with the physical ci- dimensional space, by smoothly interpolating 
the NLcrM field configuration Q into one dimension higher. The fact that 7r^_^_i{G / H) = 
Z guarantees that different ways of interpolations do not matter, and the action depends 
only on the physical field configurations in (i- dimensional space, (ii): For a. d = even- 
dimensional boundary, and for the NLaMs which describe the Anderson-delocalized 
boundary of the topological bulk, the field in the NLaM field theory is either a group 
element or a unitary matrix field (denoted hj g & G/H below - cf. last column of table 
|I]) which may, depending on the case, be subject to certain constraints. The WZW 
term responsible for the lack of Anderson localization takes on the form of a c? + 1- 
dimensional integral of the winding number density, as defined in equation (pOf ) of the 
main text below. 

To illustrate these terms, we will now give explicit examples for them in low 
dimensionalit ies : 

• d = 2: For the d = 2-dimensional topological insulator of the IQHE, the d = 1- 
dimensional NLcxM describing the edge states is the one on U(2n)/U(n) x \]{n). 
The field Q of the NLaM field theory can be parameterized a.s Q = f/^Af/ G 
U(2n) /U(n) x U(n) where U G U(2n) and A = diag (/„, — /„)• The term of relevance 
for the absence of Anderson localization at the edges of the integer quantum Hall 
insulator is oc (Jxy J dxtT[Ax] where := U'^dxU] this is a d = 1 dimensional 
analogue of the three-dimensional Chern-Simons term. This term can also be 
rewritten as a WZW type two-dimensional integral oc a^y dxdu e^^'tr [Qd^Qd^Q] 
(/i, = x,u). Here, we have extended the original d = 1 dimensional space to 
a two-dimensional region D by adding a fictitious space direction, parameterized 
by M G [0, 1]. The boundary dD of the two-dimensional region coincides with the 
original d = 1 dimensional space. Accordingly, the original NLcrM field Q{x) is 
smoothly extended to Q{x,u) such that it coincides with Q{x) when u = 0. Since 
7^2 [U(2n)/U(n) x U(n)] = Z the two-dimensional integral turns out to depend only 
on the field configuration on dD. Note that when we specialize to the case of 
n = 1, all these expressions are well-known from the coherent state path integral 
of an SU(2) quantum spin whose path integral is described by a c? = 1 dimensional 
NLaM on the target space 0(3) = U(2)/U(l) x U(l). 

• ij = 3 : For d = 3 topological insulators (superconductors) in classes AIII, Dili and 
CI, the relevant NLuMs at their d = 2-dimensional boundary (surface) are the ones 
with a group manifold as target manifold [U(n), 0(n), Sp(n) for classes AIII, Dili, 
and CI, respectively]. The relevant WZW type term can be written as a three- 
dimensional integral, J^d'^x due'^'^Hi [{g~^df^g){g~^d,,g){g~^dxg)], where g gU(?7,), 
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0(n), Sp(n) [^. Here again, the original d = 2 dimensional space is smoothly 
extended to the three-dimensional region D, in such a way that dD coincides 
with the physical d = 2 dimensional space. Accordingly, the field configuration 
g{x,y) defined on the d = 2 dimensional space is smoothly extended to the three- 
dimensional one g{x, y, u) such that g{x, y,u = 0) = g{x, y). 

• d = A : For the d = A topological insulators in classes A, All, and AI, the relevant 
d = 3 NLcrMs have target spaces G/H = G(2n)/G(n) x G(n) with G = U, O, and 
Sp respectively. The field of the NLcxM can be parameterized hj Q = WAU where 
U G G{2n). The relevant "term of topological origin" is the Chern-Simons term 
oc / d^xef'^hi [A^d^Ax + (2/3) A^A^Ax] where = Wd^U {fi = x,y,z). The 
appearance of this term at the boundary of the d = 4 topological insulators in class 
All was discussed in Ref. . 

• d = 5 : For the d = 5 topological insulators in symmetry classes AIII, BDI, CII, 
the d = 4 dimensional NLaMs have target spaces with a matrix field g which is 
an element of G/H = U(n), U(2n)/Sp(2?T,), U(2n)/0(2n), respectively. For these 
spaces the WZW term takes the form of the o? -|- 1 = 5 dimensional integral (with 
boundary) of the winding number density defined in equation (|21|) of the main text 
below. 

Note that those entries in table ^ of homotopy groups which are integers, 
7rj(G/if) = Z, have a periodicity in d equal to four. Therefore, the forms of the "terms 
of topological origin" of WZW type listed in the above low- dimensional cases will repeat, 
with the appropriate replacement of d. 



1.2.2. Terms of topological origin 0/Z2 type. The Z2 topological term for the NLaM 
at the d = 2 dimensional surface of the d = 3-dimensional Z2 topological insulators 
in symmetry class All was discussed in Refs. [^8|, Similarly, at the d = 2 

dimensional surface of the d = 3-dimensional Z2 topological insulators in symmetry 
class CII, the Z2 topological term is added to the NLcxM [pi . 



1.2.3. Classifying space. There is a third way in which the set of Cartan symmetric 
spaces listed in table |1] appear in the context of the classification of topological insulators 
(superconductors). This arises from Fermi statistics which allows for a distinction 
between those states of the (first quantized) Hamiltonian "H with energies which lie 
above, and those that lie below the Fermi energy Ep. In a topological insulator 
(superconductor) there is always an energy gap between the Fermi energy Ep and 
the eigenstates of "H lying above as well as those lying below Ep. Therefore one 
may, without closing the hulk gap, continuously deform the Hamiltonian "H into one 
in which all eigenstates below the Fermi energy Ep have the same energy (say) = — 1, 
and all eigenstates above the Fermi energy have the same energy (say) E = +1. By 
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construction, the resulting "simplified Hamiltonian" which we often denote by Q, has 
the same topological properties as the original Hamiltonian. The structure of Q for all 



symmetry classes is listed in table III (second column) of the first article of Ref. |26| for 
systems with translational symmetry. {Q will also be used again in the present article.) 
If we were to consider a "zero- dimensional" topological insulator (superconductor) where 
"space consisted of a single point" , then for each of the ten symmetry classes we obtain^ 
a matrix Q of a certain kind. It turns out that the so-obtained matrix Q for each 
symmetry class is again an element of one of the ten symmetric spaces listed in table 0, 
except that the order in which these spaces appear is different from the order in the two 
columns of table ^ The resulting assignment is listed|j]] in the last column in table [Al] . 
The symmetric spaces to which the matrices Q belong are what is called the "classifying 
space" in the classification scheme of Ref. ||2^ employing K- Theory ^. In appendix ^ 



we collect a number of interesting relationships between these three lists of the ten 
symmetric spaces: (i) the unitary time evolution operator, (ii) the NLcjM target space, 
and (iii) the classifying space. 

We end this subsection by commenting on additional information that can be read 
off from our main results, listed in Table |^. This pertains to the existence of so-called 
weak topological insulators (superconductors), as well as of zero-energy or extended 
modes localized on topological defects in topological insulators (superconductors). 

1.2.4- Weak topological insulators and superconductors. Table |^ classifies topological 
features of gapped free fermion Hamiltonians that do not depend on the presence of 
translational symmetries of a crystal lattice 0. In particular, these properties are not 
destroyed when translation symmetry is broken, e.g., by the introduction of positional 
disorder. Systems which exhibit such robust topological properties are often also referred 
to as strong topological insulators (superconductors). This is to contrast them with 
so-called weak topological insulators (superconductors) which only possess topological 
features when translational symmetry is present. As soon as translational symmetry 
is broken, such weak topological features are no longer guaranteed to exist, and the 
system is allowed to become topologically trivial. For example, systems defined on 
a (i- dimensional lattice whose momentum space is the d-dimensional torus T'^, allow 
for weak topological insulators which are not strong topological insulators. Such 
systems are topologically equivalent to parallel stacks of lower dimensional strong 

I This "simplified Hamiltonian" Q is often referred to in Ref. |^ and in the present article as a 
"projector", since it is trivially related to the projection operator P onto all filled (with E < Ep) 



eigenstates of "H by Q = 1 — 2P. (The matrix in subsection 1.2.1 that is denoted by the same letter Q 
is an entirely different object, not to be confused with Q = 1 — 2P appearing here.) 
§ From table III of the first article of Ref. [OgI, by setting the wavevector fc — 0. 



II The last row entitled "classifying space" in table Al, is nothing but a list of the matrices Q, and this 

list follows from the second column of tabic III of the first article of Ref. by letting k — 0. 

% The list of "classifying spaces" was obtained in Ref. |^ through the use of Clifford algebras. 

+ As wc will explain below, in order to classify such so-called "strong topological insulators" it is 

sufficient to only consider continuous Hamiltonians with momentum space S"^, the d-dimcnsional sphere. 
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topological insulators (superconductors). Specifically, the three-dimensional "weak" 
integer quantum Hall insulator on a lattice, discussed in Ref. is essentially a 

layered version of the two-dimensional integer quantum Hall insulator. Hence, it is 
characterized by a triplet of Chern numbers, each describing the winding of a map from 
the two-dimensional torus T^, which is a subspace of the three-dimensional momentum 
space T^, onto the complex Grassmannian, Gm,m+n(C). Similarly, in symmetry class 



AH, there exists a three-dimensional weak topological insulator [|12|, |13|, [IJ, 0, which 
consists of layered two-dimensional Z2 topological quantum states, and whose weak 
topological features are described by a triplet of Z2 invariants. 

The three-dimensional "weak" integer quantum Hall insulators, as well as the 



weak topological insulators of Refs. [12, 13, 14, 15] are both examples of c? = 3- 



dimensional weak topological insulators of "codimension" one, i.e., they can be viewed 
as one-dimensional arrays oi d = 2-dimensional strong topological insulators. The 
existence of these states can be read off from the d = 2-column and the rows of Table 
H labeled A and AH, respectively. In general, dimensional weak topological states 
can be of any "codimension" k, with < A; < d. The existence of weak topological 
insulators and superconductors in d spatial dimensions and for a given symmetry 
class can be readily inferred from table ^. That is, d-dimensional weak topological 
insulators (superconductors) of "codimension" k can occur whenever there exists a 
strong topological state in the same symmetry class, but in d—k dimensions. Specifically, 
the topological order (characterized by elements in Z or in Z2) can be specified for 
each of the (^) "orientations" of the submanifold. Moreover, it is also possible that in 
addition "strong" topological order, i.e., in the full space dimension d, exists. In the K- 
Theory description of Ref. |^ the presence of weak topological features is described by 
additional summands in the Abelian group of topological invariants (i.e., Z or Z2) which 
appears when homotopy classes of maps from the sphere (strong topological insulators) 
are replaced by maps from the torus T'^ (weak and strong topological insulators). 



1.2.5. Zero modes localized on topological defects. Another interesting application 
of table ^ is to determine whether r-dimensional topological defects in topological 
insulators (superconductors) can support localized (isolated) zero modes. For example, 
let us first consider a point-like defect (i.e., r = 0), which is embedded in a ci-dimensional 
system {d > 1) of any symmetry class. Furthermore, let us assume that at the defect 
time-reversal symmetry is broken (T = 0) but particle-hole symmetry is preserved with 
C = +1, i.e., we are dealing with symmetry class D. From table |^ we infer that in 
symmetry class Dind = r + l = l dimensions there is a Z2 classification, which implies 
that the (r = 0)-dimensional boundary can support gapless states, and, therefore, a 
point-like defect with the symmetries of class D can bind zero modes. This situation is 
realized, e.g., for a vortex core in a (p±ip)-wave superconductor, which can bind isolated 
Majorana modes. In general, whether it is possible for a r-dimensional topological defect 
of a given symmetry class to support gapless states or not is determined by the entry 
of table § at the intersection of column d = r + 1 and the row of the given symmetry 
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class. An application of this general statement can for example be found in the work 
of Ref. [Q, where a dislocation line (r = 1) in a d = 3-dimensional lattice hosting 
a weak topological insulator in symmetry class All was found to bind an extended 
gapless zero mode. (This corresponds in table |^ to the intersection of the row denoting 
symmetry class All with the column (i = r + l = l + l = 2, the latter representing the 
d = 3-dimensional weak topological insulator, which corresponds to a strong topological 
insulator in "codimension" k = 1.) To illustrate the use of table ^ we can for example 
predict that a similar binding of extended zero modes to lattice dislocation lines (r = 1) 
can also occur in weak d = 3-dimensional topological insulators of "codimension" k = 1 
in symmetry classes A, D, Dili, and C. 



1.3. Outline of the present article 

In this paper we discuss in detail the mechanism behind the dimensional periodicity 
and shift property appearing in table |^. We first demonstrate that there exists, 
for all five symmetry classes of topological insulators or superconductors, and in all 
dimensions d, a representative of the Hamiltonian in this class which has the form 
of a Dirac Hamiltonian, by constructing explicitly such a representative. We use 
these Dirac Hamiltonian representatives in the five symmetry classes with topological 
insulators (superconductors) to obtain relationships between these five symmetry classes 
in different dimensions. First we construct dimensional hierarchies ("dimensional 
ladders") relating Z topological insulators (superconductors) in different dimensions 
and symmetry classes (see sections ^ and |^). This is done by using a process of 
"dimensional reduction" in which spatial dimensions are compactified (in "Kaluza-Klein 
like" fashion), thus relating higher to lower dimensional theories. By employing the 
same dimensional reduction process, we derive in section ^ the topological classification 
of the Z2 topological insulators (superconductors) from their higher-dimensional parent 
Z topological insulators (superconductors) in the same symmetry class. Sections 0, 
H, H taken together provide a complete, independent derivation of table |^ of topological 
insulators (superconductors). 

In section |2.2| we establish a connection between the topological invariant 
(winding number) defined for topological insulators and superconductors with chiral 
symmetry and yet another topological invariant, the Chern-Simons invariant, for general 
dimensions [equation (^3|)]. For lower dimensional cases, the formula (|4^) relates the 
winding number to the electric polarization {d = 1 spatial dimensions), and to the the 
magnetoelectric polarizability [d = 3 spatial dimensions). 

In section |2]^ we list, for topological insulators in symmetry classes A and AIII the 
topological field theory describing the space time theory of linear responses in all space 
time dimensions D = d+1. For symmetry class A this is the Chern-Simons action, and 
in class AIII it is the theta term with theta angle 6 = n. (Generalizations to topological 
singlet superconductors with SU(2) spin- rotation symmetry are mentioned in section 



Finally, as already mentioned in subsection 1.2, we have discussed for all 
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dimensionalities the "terms of topological origin" of Wess-Zumino-Witten type which 
can be added to NLaM field theories on symmetric spaces. In particular, in odd 
dimensions, these are the Chern-Simons terms. 

In appendix ^ we collect a number of interesting relationships between the 
Cartan classification of generic Hamiltonians due to Altland and Zirnbauer, the NLaM 
target space of Anderson localization, and the classifying space appearing in the K- 
theory approach to topological insulators. The representation theory of the spinor 
representations of the orthogonal groups SO(A^) is briefly reviewed in appendix 0. 
Furthermore, we study in appendix |C| the influence of inversion symmetry on the 
classification of topological insulators (superconductors), combined with either time- 
reversal or charge-conjugation (particle-hole) symmetry. 

2. Dimensional hierarchy: complex case (A — AIII) 

In this section, we consider the relationship between 2n-dimensional class A topological 
insulator and {2n — l)-dimensional class AIII topological insulator. (Recall from table 
H] that Hamiltonians in symmetry classes A and AIII are not invariant under time- 
reversal and charge-conjugation symmetries (T and C), but Hamiltonians in class 
AIII possess "chiral symmetry" called S in the table, i.e., they aniz-commute with 
a unitary operator.) First we present a general discussion about class A and AIII 
topological insulators (sections p.l| and |2.2| ) and then we will focus on Dirac Hamiltonian 
representatives in these symmetry classes (sections ^.3| , |2.4| , and p^. 51 ). The insights gained 
from the study of these Dirac Hamiltonian examples are applicable to a wider class of 
insulators. Finally, in section |2.6| we discuss effective topological field theories describing 
topological response functions of class A and AIII topological insulators. 

2.1. Class A in d = 2n + 2 dimensions 

The class of problems that we will consider below is non-interacting fermionic systems 
with translation invariance. For such systems, the eigenvalue problem at each 
momentum k in the Brillouin zone (BZ) is described by 

'H{k)Mk)) = Ea{k)Mk)), a = l,...,N,,,. (6) 

Here, 'H(/c) is a A'tot x A^tot single-particle Hamiltonian in momentum space, and |u„(/c)) 
is the a-th Bloch wavefunction with energy Ea{k). We assume that there is a finite gap 
at the Fermi level, and therefore, we obtain a unique ground state by filling all states 
below the Fermi level. (In this paper, we always adjust Ea{k) in such a way that the 
Fermi level is at zero energy.) 

We assume there are (A^+) occupied (unoccupied) Bloch wavefunctions for each 
k with + N_ = Atot- We call the set of filled Bloch wavefunctions Uw^ik))}, where 
hatted indices d = 1, . . . , labels the occupied bands only. We introduce the spectral 



CONTENTS 



21 



projector onto the filled Bloch states and the "Q-matrix" by 

E 



m = EKW>(%WI. Q(k) = l-2P(k). (7) 



The ground state (filled Fermi sea) is characterized at each k by the set of 
normahzed vectors which is a member of U(A^tot), or, equivalently, by 

Q{k), up to basis transformations within occupied and unoccupied bands. Thus 
Q{k) can be viewed as an element of the complex Grassmannian Gn+,n++N-{'C) = 
U(A^+ + A^_)/[U(A^+) X U(A^_)]. For a given system, Q{k) defines a map from BZ into 
the complex Grassmannian, and hence classifying topological classes of band insulators 
is equivalent to counting how many distinct classes there are for the space of all such 
mappings. The answer to this question is given by the homotopy group vrrf(G'.m,m+n(C)), 
which is nontrivial in even dimensions d = 2n + 2. 

Define for the occupied bands, the non-Abelian Berry connection R 

A'\k)=Af{k)dk,= {u7{k)\du7{k)), /x = l,---,d, a,6=l,---,iV_, (8) 
where = —{A^^)*. The Berry curvature is defined by 

J^'^'ik) = dA""' + {A'Y' = ^Fl:t{k)dk^ A dk,. (9) 

Class A insulators in even spatial dimensions d = 2n + 2 [n = 0,1,2,...) can be 
characterized by the Chern form of the Berry connection in momentum space. The 
(n + l)-th Chern character is 

The integral of the Chern character in d = 2n + 2 dimensions is an integer, the (n + l)st 
Chern number, 

Ch„+i[J-]=/" ch„+i(^)=/ -^—^tif'^] GZ.(ll) 



'g2'*=2"+2 jg2<*=2"+2 {n + 1)! \2tt J 

Here, Jg^^ denotes the integration over dimensional /c-space. For lattice models, it 
can be taken as a Wigner-Seitz cell in the reciprocal lattice space. On the other hand, 
for continuum models, it can be taken as M.'^. Assuming that the asymptotic behavior 
of the Bloch wavefunctions approach a /c-independent value as \k\ — j- oo, the domain of 
integration can be regarded as S'^ Q 

When, d = 2 (n = 0), Chi[J^] is the TKNN integer [||, 

Chi m = ^ I tr (J^) = -L [d^k tr (F12) , (12) 



27V Jbz'*=2 27r 

* For Dirac insulators with linear dispersion, discussed below, this assumption is not entirely correct. 
If the fc-linear behavior of the Dirac spectrum persists to infinitely large momentum, the behavior of 
the wavefunctions at |/c| — > 00 is not trivial: the one-point compactification thus does not work. In this 
case, the domain of integration becomes effectively a half of 5''', and accordingly, integer topological 
numbers become half integers. However, the Dirac spectrum can be properly regularized, in such a way 
that the behavior of the wavefunctions becomes trivial at larger fc, see for example equation (|50|). This 
is also the case when the Dirac Hamiltonians are formulated on a lattice [see equation (p2|)]. 
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which is nothing but the quantized Hall conductance a^y in units of {e^ /h). The second 
Chern- number (n = 1), Ch2[-F], 

can be used to describe a topological insulator in d = 4, or, alternatively, a certain 
adiabatic "pumping process" in lower spatial dimension [|^. 



The Chern character ch„_|_i can be written in terms of its Chern-Simons form, 

Ch„+l(^) = dQ2(n+l)-lM,-^). (14) 

Here the Chern-Simons form is defined as 



Q2n+i(A :F) := - J dt tr (^J-f ) , (15) 

where 

J^t = tdA + t^A^ = tJ^ +{t^ -t)A^. (16) 

For example, 

Qi{AJ') = ^tTA, g3(A-^) = ^tr(^^d^+^^3^. (17) 
2.2. Class AIII in d = 2n + 1 dimensions 

2.2.1. Winding number. We now discuss band insulators in symmetry class Alll. By 
definition, for all class Alll band insulators, we can find a unitary matrix F which 
anticommutes with the Hamiltonians, 

{^(A;),F} = 0, F2 = 1. (18) 

It follows that the spectrum is symmetric with respect to zero energy and A^^. = A^_ =: 



A^. As a consequence of the chiral symmetry ([T8|), all class Alll Hamiltonians, as well 



as their Q-matrix (|^), can be brought into block off-diagonal form, 

in the basis in which F is diagonal. The off-diagonal component q{k) defines a map from 
BZ onto U(A^), and classifying class Alll topological insulators reduces to considering the 
homotopy group 7rd(U(A^)). The homotopy group is nontrivial in odd spatial dimensions 
d = 2n + 1, while it is trivial in even spatial dimensions, i.e., there is no non-trivial 
topological insulator in class AIII in even spatial dimensions. 

In odd spatial dimensions d = 2r2 + 1, class AIII topological insulators are 



characterized by the winding number [26 



i^2n+i[q] ■■= / U2n+l[q], (20) 

f + l 
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where the winding number density is given by 

n+l 



--l'l-l^(^) (21) 



-l)"n! / i 



CS 



2n+l 



(2r2 + 1)! \2tt 

2.2.2. Chern- Simons invariant. While derived from the Chern form, the Chern-Simons 
forms, introduced in equation (15), themselves define a characteristic class for an odd- 
dimensional manifold. This suggests that the Chern-Simons forms can be used to 
characterize also AIII topological insulators in d = 2?7, + 1 dimensions. Integrating 
the Chern-Simons form over the BZ, we introduce 

[A,T] := / Q2n+M,T). (22) 

Unlike the winding number, CS2n+i['^! is defined without assuming chiral symmetry 
and can be used for non-chiral topological insulators (superconductors). 

Before discussing how useful Chern-Simons forms are, however, observe that they 
are not gauge invariant. Neither are the integral of the Chern-Simons forms over the 
BZ. However, for two different choices of gauge A and A', which are connected by a 
gauge transformation 

A' = g-'Ag + g-'dg, r = g-'J^g, (23) 

we have 

Q2n+l{A',r) - Q2n+l{A,J^) = Q2n+l{g~^dg , 0) + d«2n, (24) 

where a2n is some 2n form [^. We note that 

n+l ^1 



Q2„+i {g-'dg, 0) = 1 (^^^ ^ dt tr {A{t' - trA'-) 

tr [((7-M(7)2"+i] f dtif-ty 
Jo 



^ , . , n+l 



= (^2n+l[g]- (25) 

This is nothing but the winding number density, and its integral 

l^2n+l[g]-= / Q2n+l{g~^dg,0) = / UJ2n+l[g] (26) 

is an integer, which counts the non-trivial winding of the map g{k) : BZ^""^^ — )■ U(A^). 
Note that 7r2n+i [U(A^)] = Z (for large enough A^). We thus conclude 

CS2„+i [A', r] - CS2n+i [A, = integer, (27) 
and hence the exponential 

W2n+i := exp{27riCS2„+i [A, T]} (28) 
is a well-defined, gauge invariant quantity, although it is not necessarily quantized. 
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The discussion so far has been generaL In particular, it is not restricted to chiral 
topological insulators. We now compute the Chern-Simons invariant for class AIII 
topological insulators in d = 2n + l. We first explicitly write down the Berry connection 
for chiral symmetric Hamiltonians. To this end, we observe that for unoccupied and 
occupied bands, the Bloch wavefucntions satisfy 

Q{k)\uf{k)) = +\nt{k)), Q{k)\u7{k)) = -\u7{k)), (29) 

respectively. Introducing 

72 




ui{k)) = ^{ 1, 6 = ±, (30) 



we rewrite (|29|) as 

q{k) 
q\k) 

We can then construct a set of eigen Bloch functions as 





(31) 



= ^ ( tm e 1 • (32) 



The A^-dimesional space spanned by the occupied states {Iw^ik))} can be obtained by 
first choosing independent orthonormal vectors which are /c- independent, and then 
from n^, 

\ul{k))n = 4^ ( , I . (33) 

From these Bloch functions, the Berry connection is computed as 

n{ul{k)\d,ul{k))^<lk, = \ {{nl\d,\n^^^e\n\q{k)\dM\k)n^^j\ dfc^ 

= \{^'re\nlq{k)\d,\q\^)n'^'\ 

= )^q{k)d,qKk%^^k,=:A% (34) 

where we have made a convenient choice (?T-a)fe — ^hh- While this looks almost like 
a pure gauge, it is not exactly so because of the factor 1/2. Note also that this 
calculation shows that the Berry connection for the occupied and unoccupied bands 
are identical. With equations ( P^ through (^^, we have succeeded in constructing 
eigen Bloch wavefuctions of the "Q- matrix", equation (0), in block-off diagonal basis, 
that are free from any singularity. I.e., we have explicitly demonstrated that there is no 
obstruction to constructing eigen wavefuctions globally. We emphasize that this applies 
to all symmetry classes with chiral symmetry (AIII, BDI, Dili, CII, CI) for any spatial 
dimension {d = 2n as well as d = 2n + 1). 

Alternatively, one can construct another set of eigen Bloch functions 

i„,.),._L(^.«"iy ,35, 
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For these Bloch functions, the Berry connection is computed as 

s{ul{k)\d,ul{k))sdk, = \ {^{nlq\k)\dM{^)n'^ + Kl^.l^^)] 

= \{^knlq\k)\dM{^)n'^^^\ 

= \\^\k)d,q{k)\.^k,=:A% (36) 

where we have made a convenient choice = ^-^i- The two gauges are related to 

each other by 

= g-W'g + g~'dg, (37) 
where the transition function g is the block off-diagonal projector, g = q: 

q\k) 



^q{k)d,q\k) q{k) + q\k)d,q{k) = ^ [d,q\k)\ q{k) + q\k)d^q{k) 



2 

= ~lq'{k)d,q{k) + q^{k)d,q{k) 

= lqKk)dAk)- (38) 

We now compute the Chern-Simons invariant for class AIII topological insulators 
in c? = 2n + 1. In the gauge where A = A^ = (l/2)(g~-'^dg), we have 

d^= id(g"Mg) = -i(g-Mg)(g-Mg), 

=tdA + t'A' = (-^- + pi (g-Mg)2. (39) 



Then, 



Q2n+i{A, ^) = (^^^ ^ dt tr {AJ^n 

1 r i V' f f t evdt 



, , tr \{q-^dqY''+^ 
n\\2Ti /n V 2 4/2 '^^ J 



= ^ (^) ^ - l)"d« tr [(g-Mg)2-+^] , (40) 

which is a half of the winding number density, 

Q2n+M^J^^) = ^U2n+l[q]. (41) 

Hence we conclude 

CS2„+i [A^J^""] = f Q2n+l{A^J^'') = lu2n+l[q]. (42) 

As a corollary, 

1^2^+1 = exp{27ri CS2n+i [A, J^]} = exp{7ri z/2„+i[g]}. (43) 

Thus, while in general the quantity W2n+i is not quantized, for class AIII Hamiltonians 
in d = 2n + 1, W2n+i can take only two values, W2n+i = ±1- 
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Physically, CS2n+i [-^7 -F] in = 1 (n = 0) spatial dimension takes the form of 
the U(l) Wilson loop defined for BZ"'^^ ^5^. It is quantized for chiral symmetric 



systems Also, the logarithm of Wi represents the electric polarization |57, 58 

For a lattice system, the non-invariance of CSi J-"] [i.e., it can change under a gauge 
transformation, CSi [^,-7-"] — )■ CSi [^, J-"] + (integer)] has a clear meaning: since the 
system is periodic, the displacement of electron coordinates has a meaning only within 
a unit cell, and two electron coordinates that differ by an integer multiple of the lattice 
constant should be identified. 

In d = 3 (n = 1) spatial dimensions CS3 represents the magnetoelectric 
polarizability [pSj |59| , |60| . While this was discussed originally for three-dimensional 
Z2 topological insulators in symplectic symmetry (class All), the magnetoelectric 
polarizability is also quantized for class AIII. 

2.3. Dirac insulators and dimensional reduction 

We start from a.d = 2n+3-dimensional gapless Dirac Hamiltonian defined in momentum 
space, 

d=2n+3 

^ST(^)= E ^'^r^^n+s)- (44) 

a=l 

Here, fca=i,...,d are (i- dimensional momenta, and r"2n+3)^"'''^ (2"+^ x 2"+^)-dimensional 
Hermitian matrices that satisfy {^l2n+3)^^\2n+3)} ~ '^^a,b- Some results for the gamma 
matrices, on which the following discussion is based, are summarized in appendix p. 
The massless Dirac Hamiltonian 'H'^^^^^^-^ {k) cannot be realized on a lattice in a naive 
way, because of the fermion doubling problem. 

By replacing k2n+z by a mass term, we obtain & d = 2n + 2 dimensional class A 
topological Dirac insulator, 

d=2n+2 
a=l 

The topological character of this Dirac insulator will be further discussed below. 

By setting in addition /c2n+2 = 0, we obtain an insulator in one dimension lower, 

d=2n+l 

^&T(^'H= E ^«r^2„+3)+^rj;+3)- (46) 

a=l 

By construction, ^^l^n+t)^ ^ ^) anticommutes with T^^^g^, 

'Htn:t)\k.ml^^S,)]=^. (47) 

and hence it is a member of class AIII. 

This construction of the lower-dimensional models from their higher-dimensional 
"parent" is an example of the Kaluza-Klein dimensional reduction. To obtain 
^(2n4^3)^(^' ''^) froni 'H^2^^_|^3)^(^' ''^)' ^^^^ compactify the (2n + 2)-th spatial 

direction into a circle S^. The (2?7,+2)-th component of the momentum is then quantized. 
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^2n+2 = 27rAi'2n+2/^, wheie A^2n+2 ^ ^5 and £ is the radius of . The energy eigenvalues 
now carry the integral label N2n+2, in addition to the continuous label, A;j=i^..._2n+i- By 
making the radius of the circle very small, all levels with N2n+2 7^ have very large 
energies while the levels with A^2n+2 = are not affected by the small radius limit. In 
this way, in the limit £ — )■ 0, there is a separation of energy scales, and we only keep 
the states with N2n+2 = 0, while neglecting all states with N2n+2 7^ (the Kaluza-Klein 
modes). 

Similarly, 'H'l^2n+3)^(^^^) naturally obtained from 'H'^2n+3)^ dimensional 
reduction, if we also include a U(l) gauge field. We can introduce a (external, 
electromagnetic) U(l) gauge field aj=i,...^2n+3 by minimal coupling, 'H'(2n+3)^ ~^ 



'-(2n+3) 

^(2rf+3)^(^ + a). In the compactification of the (2n + 3)-th spatial coordinate, we 
again keep modes with k2n+3 = only, while neglecting all the Kaluza-Klein modes with 
^2n+3 7^ 0. The vector field (gauge field) aj=i^...2n+3 in d = 2n + 3 dimensions can be 
decomposed into a vector field (ai=i,...2n+2) and a scalar field (02^+3 = (fi) in d = 2n + 2, 
leading to 

d=2n+2 

+ ^) = E + «^)r(2n+3) + (48) 

i=l 

Switching off the d = 2n + 2-dimensional gauge field, and assuming the scalar field is 
constant cp = m, we obtain '^'^(^^^^■^ {k , m) . 

For the massive Dirac insulator (^5)) in class A in 2n + 2 dimensions, the n-th Chern 
number Ch„+i is nonvanishing [^. Accordingly, for the massive Dirac insulator (^) 
in class AIII in 2n + 1 dimensions, the winding number and the Chern-Simons form are 
non-zero. 

As the Chern-Simons form can be derived from its higher-dimensional parent, the 
Chern form, one may wonder if there is a connection between CS2n+i = |^2n+i of a class 
AIII Dirac topological insulator in d = 2n -|- 1 and Ch„_|_i of a class A Dirac topological 
insulator in d = 2?7, -|- 2. We now try to answer this question. 

Let us start from 

d=2n+2 

^&T(^'H= E ^'^r^2n+3)+MA;)r?2':S3). (49) 
a=l 

Here we have regularized the Dirac spectrum by making the mass term k dependent, 
m{k) = m — Ck^, with sgn(C) = sgn(m). (50) 

With this regularization, as \k\ — )• 00, 'H'l^2n+3)^(^^''^) ^ -'^(2^+3)' ^^'^ hence the 
wavefunctions become fc- independent (i.e., the BZ is topologically 5"*^). The sign choice 
of the constant C, sgn(C) = sgn(m), makes the Chern invariant for the Dirac insulator 
non-zero. 

By dropping k2n+2, we get the (regularized) Dirac insulator in one dimension lower, 

d=2n+l 

-^tn^t'ik.m) = E kaTt2n+,) + m{k^)T^^Zy (51) 
a=l 
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where k = (/ci, . . . , /c2n+i)- As mentioned before, 'H'^^^^^^^{k,m) is chiral symmetric, 
{^?2;?3y(^'"^)'r(2"J"+3)} = 0. We work in the basis where r^^';;^^^^^ is diagonal In this 
basis, block-off diagonal, 

«ffS'(*.™)=(^,(^ T). (52) 

Correspondingly, 

njd=2n+2(,-. _ ( D{k) \ 

where A(/c2„,+2) = /^2n+2, and D{k) = D{k,k2n+2) satisfies 

D{k,k2n+2) = D{k), when /C2n+2 = 0. (54) 

We now look for eigenfunctions of 'Hf^n+t'^ with negative eigenvalues, 

-.XL) -<^-' ( : ) ■ 



To this end, we first solve the following auxiliary eigenvalue problem 

= -m ( M , (56) 



D{k) 
D\k) 





where A = \J k"^ + [m{k)]'^. We see from ( p3D and ( pSj) that the solutions are given by 

V2[qKk)n,) V2[ -n, )^ ^^^^ 

where q{k) and q^{k) are the off-diagonal blocks of the Q-matrix derived from the 
auxiliary Hamiltonian. If we specialize to the case where k = {k,0), the projector 
'^^ '^(2rf+3)^(^) given in terms of q{k) and q^{k). We then construct two different sets 
of normalized eigenfunctions of "H^g^^s)^ 



and 




K(^))N^^7__i 1, (58) 



\u7{k))s= , ^ I (^T^)^-^ ), (59) 

" v/2A(A- A) I Ana ' 



with the eigenvalue —A = — a/A^ + A^. When specialized to /c = {k,0), these 
wavefunctions yield the eigenfunctions (0) of the chiral Dirac Hamiltonian H'f^n+t)^ ■ 

Since the Hamiltonian 'H'^^^"^^(fc) is characterized by the non-zero Chern number, 
Ch„+i [J^] 7^ 0, it is not possible to define wavefunctions globally We thus need to split 
the BZ (BZ'^=2"+2) into two parts, BZ^=2"+2 and BZ^=2"+l We choose the interface 
of these two patches as 9BZ^=2"+2 = -9BZ^=2"+2 ^ ^r^d=2n+i^ ^j^j^j^ ^ ^^^^^ 
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Figure 2. (a) The d = 2n + 2 dimensional Brillouin zone, BZ'^-2"+^ - S'2"+2 = 
^2«+i ^ g^^^^ its d = 2n + 1 dimensional descendant BZ'^=2"+^ - 3'^'^+^ located at 
the equator of S'^"'*'^ . (b) Splitting the c? = 2n + 2 dimensional Brillouin zone in half, 
or embedding d = 2n + 1 dimensional Brillouin zone into a higher-dimensional one. 



dimensional Hamiltonian "lives". For each patch, BZ^"g^"''~^, we can choose a global 
gauge. Observe that the wavefunction |m^(/c))n is well-defined for /c2n+2 > 0, and has a 
singularity at (fc, fc2n+2) = (0, —\frnfC\ Similarly, |M^(fc))s is well-defined for ^2^+2 < 
and singular at (A;, fc2n+2) = (0, \frnfC). 

The two gauges are thus complementary and at the boundary c^BZ^^^"'*'^ = 
— 9BZg^^"^^ they are glued by the transition function 

= g-'A^'g + g~'dg, (60) 

where g{k) G U(iV_) with k G (9BZ^=2'"+l Then, 

Ch,+l[^]=/ Ch„+i(^)= / dQ2n+l(A-F) 

dg2n+l(^'', ^) + / dg2n+l(^', ^) 

J=2n+2 I T:,r,d—2n+2 

= / Q2n+liA'',J')- f Q^n+liA^J'), (61) 

where we have used Stokes theorem. From formula (0), it follows that 

Ch„+i[J-]=/ Q2n+i{9''dg,0). (62) 

J9Bz^r^"+^ 



By use of equation ( |25D we conclude that the RHS of the above equation is nothing but 
the winding number. Thus, the Chern number can be expressed as a winding number 
of the transition function g. For the case we are interested in, i.e., when a = 2n + 1 
topological insulator "lives" on c?BZ^^^"^^, the transition function is given by the off- 
diagonal block of the projector, g{k) = q{k). This is how two topological insulators in 
d = 2n + 2 and d = 2n + l are related. 

2.4. Example: d = "i ^ 2 ^ I 

Consider three mutually anticommuting, Hermitian matrices 

pa=l,2,3 ^ 1^^^ J^g2^| 
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where (Jx,y,z are the 2x2 Pauh matrices. These matrices can be used to construct a 
d = 3 gapless (Weyl) chiral fermion 

3 

^(3f(^) = E^«r^3)- (64) 

a=l 

By replacing the third component of the momemtum by a mass term, ka=z — )■ m, we 
get a. d = 2 dimensional Dirac Hamiltonian, 

2 

a=l 

= k^a^ + kyay + ma^. (65) 
This is a class A Hamiltonian. The Bloch wavefunctions are given by [A(A;) := 




(66) 

We note that \u~{k)) is well-defined for all k when m > 0. (When m < 0, by choosing 
a gauge properly, we can obtain a similar well-defined wavefunction.) Assuming m > 0, 
we adopt |m~(/c)), for which we obtain the Berry connection 

= +2A(fT^- '^»'^-'"' = -2A(fT;^)- («^' 

and the Berry curvature 

vm 

E;,y {k, m) = dk^Ay - 9^v4^ = - — . (68) 
The Chern number is non-zero. 



au|^l = -yd^«:F., = -yd'»:^ = 5^. (69) 

which is nothing but the Hall conductance a^y This model can be realized, at low 



energies, in the honeycomb lattice model introduced by Haldane p2 



Finally, setting /c2 = 0, we get a d = 1 dimensional Dirac Hamiltonian, 

T^fgf (/c, m) = k^a^ + ma^. (70) 

This is nothing but the chiral topological Dirac insulator in class AIII, as it anticommutes 
with ay. A topological invariant can be defined for the d = 1 dimensional Dirac 
Hamiltonian as follows. We first go to the canonical form by {a^, ay, a^) — >■ {a^, cr^, — o'y), 
^(3)^(^5 m) — )■ kxax — may. This d = 1 dimensional Hamiltonian describes the physics of 
polyacetylene (see, for example, Refs. [^, ^If-) The Bloch wavefunction with negative 
eigenvalue is 

V2(*i + m?) \ \ I 
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The off-diagonal block of the projection operator Q{k) in this basis is 

q(^k) = - ^""^'^^ (72) 

and the winding number is given by 

i /" 1 , i f°° —im 1 m , 

2tt Jbz 27r J^,,^ k^ + m^ 2 |m| 

From equation ( [TTD the Berry connection is obtained, 

^(A;,) = {u~ik,)\du-{k,)) = -——^dk,. (74) 

Then, we find 

This is a half of the winding number z/i[g] as expected. Hence, the Wilson "loop" is 
given by 



Wi = exp {27riCSi [A, J^]} = exp / tr ^ 

Jbz 



,±7ri/2 



(76) 



Here, we mention, however, that there is a subtlety in computing the Wilson loop. 
In the basis where the Hamiltonian in momentum space is real, 'H'^^^{k, m) = kxfJx+fnaz 
([70|), the Bloch wavefunction is real for a given /c, and as a consequence the Berry 
connection A{kr^) vanishes identically. One would then conclude Wi = 1, not Wi = — 1, 
although ui = |sgn(m). This puzzle can be solved by properly regularizing the Dirac 
insulator. For simplicity, let us replace k^ by sinfc^., and m by m — 1 + cos{kx). With 
such a regularization, the BZ is topologically S^, and one finds a singularity in the 
wavefunction at kr^ = ir, where the phase of the wavefunction jumps by it (i.e., there 
is a Dirac string at k^ = vr.). If, on the other hand, we use a different basis (|7TD , we 
can avoid to have a Dirac string in the BZ, and the Wilson loop can be obtained by 
integrating the Berry connection over the BZ. 

2. 5. Example: = 5 4 -> 3 

Consider five mutually anticommuting, Hermitian matrices 

pa=l,...,5 ^ 1^^^ -i/37^}, (77) 

where we are using the Dirac representation, 

-(:o)' "-['^-^ ^^^(-)- 

Observe that axOiyazP = — i/375. These matrices can be used to construct a. {d = 5)- 
dimensional gapless chiral fermion 

5 

^(5fW = E^«r^5)- (79) 

a=l 




Figure 3. (a) Winding number 1^3 for Hamiltonian ( |8^ ) as a function of 7715. (b) Two- 
dimensional energy spectrum of the surface states of model ( |8^ ) with mass TO5 = +0.5. 
There are two inequivalent surface modes in agreement with the winding number 
i^3{m5 = +0.5) = -2. 



By replacing the fifth component of the momemtum by a mass term, ka=5 — )■ m^, we 
obtain a. d = 4 dimensional Dirac Hamiltonian, 

4 

nl^^'ik, ms) = J2 f'-^h + ^5r?5). (80) 

a=l 

It is known that for this gapped Hamiltonian, the second Chern number is non-zero 

Finally, setting = 0, we get a d = 3 dimensional Dirac Hamiltonian, 

3 3 
Hf^fik, ms) = ^-^h + ^5^h = E - ^5i/37'- (81) 

a=l a=l 

This is nothing but the chiral topological Dirac insulator in class AIII discussed in Ref. 



26| . The tight-binding version of this model for a simple cubic lattice is given by 

nfffik, ms) = ^-^h +(^5 + Yl ^-j ^h- (82) 

a=l \ a=l J 

By taking open boundary conditions in the z direction and periodic boundary conditions 
in the x and y directions, we can study the surface states of this model. When the 
winding number z/3 is non- vanishing, there are l^sl surface Dirac states, which cross the 
bulk band gap (see figure 

We now study the Berry connection of the above class A {d = A) and class AIII 
{d = 3) Dirac insulators. To treat these two cases in a unified fashion, let us consider 
the following d = 4 Dirac Hamiltonian 

3 

Tiik, mo) = Y + "^0/3 + kw{-if3'y^). (83) 

a=l 
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Upon identification ky^ — it can be also viewed as a c? = 3 dimensional Dirac 
Hamiltonian with two masses. The four eigenvalues are 

E{k) = ± A(A;), A(A;) = ^Jk^ + mg. 
The two normalized negative energy eigenstates with E{k) 

1 



4) 



-\{k) are 



\u^{k)) 



v/2A(A + mo) 



iky, -\- kz 







y A + mo / 



\u^{k)) 



v/2A(A + mo) 



v 



A + m^o 




^5) 



The two normalized positive energy eigenstates with E{}i) = +X{k) are 

1 



v/2A(A - mo) 



-ik,,,, ky 







\ut{k)) 



\ A - mo J 

kx + iA^y 



v/2A(A - mo) 



v 



A — mo 




(86) 



Note that if mo > 0, |^ii'2(^)) are not well-defined at \{k) = mo (i.e, k = 0). Observe 
that when mo > 0, \ui2{k)) have the form of wavefunctions expected from the general 
consideration (corresponding to the gauge choice). 

When mo = 0, the Hamiltonian anticommutes with f3 and is block off-diagonal. 
The off-diagonal component of the Q-matrix is given by 



qik) = ^{k- a - ims) 



(87) 



By noting (yU = 1, 2, 3) 



and also, 



k 1 . —1 

9^(1 = ^{k- a- ims) - -^c^M' = ^ i^t, - k ■ aa^ - ima^) , 



dU-e^^^^r [q^d^qq^d^qq^dpq] 



d?k tr [-3im5(A; ■ a)a^{k ■ a)au(Tp + imlapa^ap] 



2i 



-3im5)e^^p— / d k 



12m, / d^k 



d^k 



(F + m2)3 



(A;2 + m2)2' 



(89) 
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the winding number is given by 

12 



247r2 



dk 



47rfc2 



(P + m2)2 



1 vr mg 

r47r X 



1 ms 



27r2 



4 Imsl 2 \m^\ 



(90) 



We now compute the Chern-Simons invariant (magnetoelectric polarizabihty) for 
the massive Dirac Hamiltonian. For the lower two occupied bands, we can introduce a 

{ul{k)\dul{k)) {a,h = 1,2). The U(2) gauge field is 



U(2) gauge field by A1f{k)dk^ = 
given by, in the matrix notation 

i 



A. 



A, 



A. 




2A(A + mo 



2A(A + mo) 



kw ~l" i^2 
k,. 



' krf 



i/c,. 



ikyj kz -\~kx 
kw '^kx ~t~ ky 



2A(A + mo) \ ikx + ky 



kn 



(91) 



The gauge field can be decomposed into U(l) {or) and SU(2) (a-' '^'^'^) parts as 



AM = a%kf-^ + aUk)^. 



The U(l) part of the Berry connection is trivial, whereas the SU(2) part is given by 



(92) 



ko. 



-kz 



~\~ky 



A(A 






'kz 


A(A 


+ mo) 


-i 


-ky 



a. 



y 



A(A + mo) ' 
— h 



1- 



y 'A(A + mo)' 

1 



'A(A + mo)' ^A(A + mo)' 

The Chern-Simons form can be computed as 

-1 f „.,„ / . „ . 2 



A(A + mo) ' 

A(A + mo) ' 
—h 

A(A + mo) ' 



(93) 



CS, 



«7r^ 

1 
8^ 

TX Imsl 



d=^A;e^^^tr ( A^d.A, + ^A^A.A, 
d'k- 



3 m5(2A + mo) 



A3(A + mo)^ 



■arctan 



[mg + mo 



_ 



mo 



{ml + miy/'^ + mo 



1/2 



When mo = 0, 



CS, 



1 1715 



4 jmsl 



(94) 



(95) 



As the CS term is determined only modulo 1, the fractional part |sgn(m5) is an intrinsic 
property. On a lattice, as there are two-copies of the 4x4 Dirac Hamiltonians, the total 
Chern-Simons invariant is given by CS3 = ^sgn(m5), and e^'^^^^^ = —1. 

There is a subtlety in computing CS3 and W3, which is similar to the case of the 
Wilson loop Wi- When mo 7^ and = 0, the Chern-Simons form is zero identically, 
and one would conclude W3 = 1, while chiral symmetry for the case of mo 7^ and 
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ms = allows us to define 1/3, and z/3 = |sgn(m). As before, this puzzle can be solved 
by properly regularizing the Dirac insulator. With such a regularization, the BZ is 
topologically S"^, and one finds a singularity in the wavefunctions. If, on the other hand, 
we use a different basis, we can avoid to have a Dirac string in the BZ. 



2.6. Description in terms of field theory for the linear responses in D = d + 1 
space-time dimensions 

For class A and AIII topological insulators or superconductors, there is always a 
conserved U(l) quantity; either electric charge or one component (the ^-component, 
say) of spin is conserved |6^. Transport properties of these conserved quantities in 



these topological insulators or superconductors in d spatial dimensions can be described 
in terms of an effective field theory for the linear responses in D = d + 1 space- 
time dimensions, which can be obtained by coupling an external space-time dependent 
gauge field to the topological insulator or superconductor (we will choose the Dirac 
representative; see the discussion near (^)), and by integrating out the gapped fermions. 
One can then read off the theory of the linear responses from the so-obtained (classical) 
effective action for a^. 

The theory for the linear responses of a topological insulator in symmetry class A 
and in d spatial dimensions is then given by the Chern-Simons action in D = d + 1 
space-time dimensions, 

'^cs='"'"'H = 2vriA; / Qsn+iH, k = Ch„. (96) 

J^D = 2n + l 

Here, Q[a] is the Chern-Simons form of the U(l) gauge field defined in D = 2n + 1 
dimensional (Euclidean) spacetime, and the level k of the Chern-Simons term is given 
in terms of Ch„, which is the Chern invariant computed from the Bloch wavefunctions 
in d = 2n space dimensions. 

On the other hand, the theory of linear responses of a topological insulator in 
symmetry class AIII in d = 2n — 1 spatial dimensions is given by a theta term in 
D = 2n space-time dimensions. 



Sl/='^[a] = iO / ch„[/], 6 = 7rz/2n-i, (97) 

where / is the field strength of the external gauge field a^, and the theta angle 9 is related 
to the winding number V2n~ii which can be computed from the Bloch wavefunctions. The 
integration of fermions in D = 2n dimensions is trickier than in D = 2n + l dimensions; 
while a simple derivative expansion of the fermion determinant would appear to yield 
the Chern-Simons action (p6|), a more careful evaluation yields in fact the theta term. 



which is obtained from the Fujikawa Jacobian pH . 

The theta angle in the theta term (|97D can take on a priori any value, if there is no 
discrete symmetry that pins its value. In class All Z2 topological insulators in D = An 
dimensions, time-reversal symmetry demands 9 to be an integer multiple of vr. Similarly, 
in class AIII, chiral symmetry pins 9 to be an integer multiple of tt. This can be seen 
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from the symmetry properties of the Chern-Simons and theta terms, 

'y- . gD=2n+l ^ ^ j^-jn^D=2n+l gD=2n y ^ '^^n+l^D=2n 



n . QD=2n+l 



CS ^ ^ti ' \ ^9 ) 

n+lQD=2n+l cD=2n , I T\noD=2n 
1 "Jfl 



^ (-l)"+^5^s='"+\ ^e^='" ^ (-1)"5,^=^". (98) 

Here, the (anti-unitary) time-reversal symmetry operation T may either square to plus 
or minus the identity. Similarly, the charge-conjugation symmetry operation C may 
either square to plus or minus the identity. Since the chiral symmetry S can be expressed 
as the product of these two symmetry operations, i.e., S = T ■ C, one obtains 

S : S^=^'' -5e^=2". (99) 

Together with the 27r periodicity of 6, it follows that 6 is either or tt. 



3. Dimensional hierarchy for Z topological insulators: real case 

3.1. Dimensional reduction 

We now consider eight "real" symmetry classes. In particular, we consider topological 
insulators (superconductors) with integer classification, Z or 2Z. (Z2 cases will be 
discussed in section ^.) A connection similar to the one between classes A and Alll can 
be established. 

The massive Dirac Hamiltonians relevant to this section are all obtained from the 
following massless Dirac Hamiltonian 

d=2n+3 

(0): ^?2;+T(^)= E ^«r^2n+3), (100) 

a=l 

by replacing some momenta by a mass, or by simply dropping them. Specifically, we 
will consider 

d=2n+2 



[III] 



^(2n+3) 



^'{2n+t)^^^^^ 



{2n+3) '''■J- (2n+3) ' 



{2n+3) ^ '"J- (2n+3)' 



{iv) : 



a=l 
d=2n+l 

a=l 
l=2n 

^(2n+3)(^, m)=Yl ^'^r^2n+3) + ^M, 
a=l 

d=2n-l 

^&V(^>H= E ^ar^2n+3) + ^-^ 



a=l 



where 



\y( ._ -^2^+3 p2n+2 p2n+l 
JVl .— li (2„+3)J- (2n+3)^ (2n+3)- 



(101) 
(102) 
(103) 
(104) 

(105) 



The massive Dirac Hamiltonian ^^2rf+3)^(^' obtained from its gapless parent 
7^^2^^3)^(^) by replacing /c2n+3 with a mass m. By further removing k2n+2, we obtain 

d=2n+1 
(2n+3) 



'HfZ^?t\k,m). This is essentially the same procedure as in the previous section 
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Table 4. 



Symmetry class of Dirac Hamiltonians '^(^2n+3)^(^^ (100), 



n 



(2n+3) I'^i'^J 



(IMl), np^^'{k,m) (|0|), Hf2f+"3)(fc,TO) (|10|) and, nl-^^^-^^k , m) marked 
by "0" , "i" , "m" , "in" , and "iw" , respectively. Chiral and non-chiral symmetric classes 
are colored in blue and red, respectively. 



where we obtained the class AIII topological Dirac Hamiltonian from its parent class 
A Hamiltonian. On the other hand, ^^2n4*3)(^' '^^^ 

be obtained from H^2n?3)^(^) by 
first removing three components of the momentum k2n+3,2n+2,2n+i, and then adding Ai 
as a mass term. Finally, n^^^-^^k, m) is obtained from 'H^2ri4^3)(^' ''^) dimensional 
reduction. 

Here, an important difference from the dimensional reduction in the complex case 
is the shifting of symmetry classes. The parent Hamiltonian Hl^n+t)" (and 
only one) of four discrete symmetries, T (squaring to plus or minus the identity), or C 
(squaring to plus or minus the identity), depending on the dimensionality d (see table 
§ and appendix |B|). That is, it is a member of any one of AI, D, AH, and C. 

3.1.1. From d = 2n + 3 to d = 2n + 2. When dimensionally reducing 'Hpn+s)^ (^) 
^f2rf+3)^(^' '^)' symmetry of the parent Hamiltonian 'H'^^^n+s)^ broken by the 
mass term: While the sign of the kinetic term is reversed under the discrete symmetry, 
the sign of the mass term remains unchanged. Thus, the new massive Hamiltonian is a 
member of a different symmetry class. By this procedure, we get an even dimensional 
Hamiltonian 'H^2^^^3)^(^' ''^) with the "shifted" symmetry, 

AI^C, C^AH, AH^D, D ^ AI. (106) 

This can be easily proved by taking an explict form of gamma matrices (see appendix 
Observe that, while we have chosen to remove /c2n+3 and hence T^^^g^, any other 
component could have been removed instead of /c2n+3 and T^^^g^. However, all these 
different choices are unitarily equivalent, as they are simply related by a permutation 
of Clifford generators, r^2„_,_3). 

3.1.2. From d = 2n + 2 to d = 2n + l. Further dimensionally reducing ^('2rf+3)^(^' ^) 
^(2n4S)^(^' obtain a chiral symmetric Dirac Hamiltonian. Together with existing 
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discrete symmetry, this chiral symmetry yields one more discrete symmetry. One can 
see the symmetry of the newly generated Hamiltonian is shifted, 

D ^ BDI, All ^ DIII, C ^ CII, AI ^ CI. (107) 

This again can be easily proved by taking an explict form of gamma matrices (see 
appendix P). 

3.1.3. From c? = 2n + 3 to d = 2n. Let us now consider 'H^^^3)(fc, m). One can check 
that this Hamiltonian belongs to the same symmetry class as its parent 'H'l^2n+3)^ (^) ■ 
To study the topological properties of 'H('^^_[^3-)(/c, m), observe that, by construction, 
'^(2rf+3)(^' '"^) commutes with a product made of any two of 'r'^(^2n+3)^~^'^''^^~^'^ i 

[ntnUk, m), r^2n+3)r?2n+3)] = (108) 

where (a, b) = {2n + 1, 2n + 2), (2n + 2, 2n + 3), or (2n + 3, 2?i + 1). This suggests that 
the Hamiltonian can be made block- diagonal. We note that in the representation of the 
gamma matrices we are using, the Hamiltonian is block diagonal, 

d=2n 
a=l 

By use of a unitary transformation for the lower right block, 

r?2J;\) rj2„+i) = ^(2=n+i)(^, -m), (110) 

'^{2n+3)(k,iTT') reduces to two copies of H'l^2n+i)^k, —m). On the other hand, we have 

already seen that each 'H^2n+i)(k, —rn) = J^tzT kaJ'^2n+i) ~ ''^^f2n+i) ^ non-trivial 
Chern number ±1. We thus conclude that the Chern number characterizing Hf2n+3)^k) 
is ±2. 



3.1.4. From d = 2n to d = 2n — 1. Finally, we discuss the topological insulators 
(superconductors) 'H'l^2n+3)^(k,m) marked by "-iu" in table ^. One can check that the 
symmetry class of this Hamiltonian is shifted from that of '^f2n+3)^k,m) as 

AI^CI, C^CH, AH ^ DIII, D^BDI. (Ill) 

Just as 'H'^2n+3)iki^) is characterized by a non-zero Chern invariant which is even, the 
descendant 'H'l^2n+3)^(k,fn) is also characterized by an even winding number. 

The symmetry properties of all the Dirac Hamiltonians discussed here are 
summarized in table ^. 

3.2. Example: d = 3^2^1 

Consider a. d = 3-dimensional gapless Dirac fermion 

nik) = k^a^ + kyay + k^a^. (112) 
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This Hamiltonian is a member of class All, since 

ay'H{k)ay = H*{-k). (113) 

Starting from this Hamiltonian, we dimensionally reduce this Hamiltonian successively 
and obtain the lower dimensional Hamiltonians with shifted symmetry classes, AH — >■ 
D — BDl. We replace kz by a mass term, and consider a. d = 2-dimensional massive 
Dirac fermion 

Tii^k, m) = kxdx + kyCy + mcTj. (114) 

This Hamiltonian is now a member of class D since 

axT-L* {k,m)ax = —'H{—k,m). (115) 

As mentioned before, the Chern number is non-zero for this model, Ch.2 = |sgn(m). 
This Hamiltonian is in the same universality class as the chiral p + ip-wave topological 
superconductor [|, E%. 



The single-particle Hamiltonian ( |114|) , of course, looks identical to (|65|). This 



degeneracy is lifted once we consider more generic types of perturbations than the 
mass term, which can in principle be spatially inhomogeneous. Also, because of the C 
symmetry [see equation ( |115| )], the Hamiltonian ( |1 14| ) can be viewed as a single-particle 



Hamiltonian for real (Majorana) fermions, while this is not the case for a class AIII 
single-particle Hamiltonian. In other words, a class D system can be written in the 
second quantized form as 

I E {ala_,)nk,m)h>^ Y (116) 



H 

2 

fceBZ- 



2 \ -~-k 



where Ci\/a^ is a fermion creation/annihilation operator with momentum k. 

By further switching off ky, or by the dimensional reduction in the ^/-direction, 

T-L^kym) = kxCTx + maz. (117) 
This Hamiltonian is a member of class BDl since 

axTi* {k,m)ax = —'H{—k,m), azT-i*{k^m)az = l-L{—k,m). (118) 
This model can be realized as a continuum limit of the lattice Majorana fermion model 



discussed in Ref. pS|. For this d = \ topological massive Dirac Hamiltonian, the 
winding number and the Wilson loop can be computed in the same way as we did for 
the d = 1 class AIII topological Dirac insulators (|TD|). From ([75D, we immediately see 
vi = Isgn(m). 

3. 3. Example: d = 5 ^2 

Consider = 2 ■ 1 + 3 dimensions. We can take the following five matrices as gamma 
matrices, 

pa=l,...,5 ^ ^ ^^^^ ^ ^ ^^-^-^g^| 
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Note that 

^y{^hY^y = +n^y (120) 
Thus, the gapless Hamiltonian 

5 

n{k) = Y^kaTl,^ (121) 

a=l 

is a member of class C, as seen from 

a,A;,(r^5/a, = -(-Ur^5)- (122) 

We now remove three momenta, and add one mass term. There are four possible mass 
terms that are compatible with the class C symmetry 

Tx.zfi ® Ty ® CTq. (123) 

To construct a kinetic term we pick two gamma matrices from the list ( |119| ). For 
any choice of two gamma matrices there is only one mass term in equation ( |123| ) that 
anticommutes with the chosen kinetic term. In this way, we obtain, for example, the 
following massive Hamiltonian 

l-iik^ m) = Ty® a^kx + Ty® a^kz + rriTQ ® ay. (124) 

Observe that i{Ty® ay){Tx® o'q){tz® (Jo) = TQ®ay. It is easy to see that this Hamiltonian 
has the doubled Chern number, Chi = ±2. To see this, first rotate Ty Tz, 

l-iik, m) ^ Tz® CTxkx + Tz® (Jzkz + rriTo ® ay 

CTxkx + <Jzkz + CTym 

—(Txkx — azkz + ayUi 

and then apply a unitary transformation. 



(125) 





(^xkx 



axkx + crzkz + (Jym 

The d = 2 dimensional Hamiltonian constructed here has the same topological 
properties as the d = 2 dimensional d + id- wave superconductor discussed in Ref. [^]. 



(126) 



4. Z2 topological insulators and dimensional reduction 

In the previous sections we have studied Z topological insulators (superconductors), 
which are characterized by an integer - the Chern or winding number. In this section 
we discuss Z2 topological insulators (superconductors) and show how they can be derived 
as lower dimensional descendants of parent Z topological insulators (superconductors) 
(see table |]). Recently, Qi et al. |2^ have employed this procedure to derive descendant 
Z2 topological insulators (superconductors) for those symmetry classes that break chiral 
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symmetry. Here, we show how this approach can be generahzed and apphed to all eight 
symmetry classes of table ^ We first treat chiral symmetric systems, and then review 
the case of chiral symmetry breaking classes. The reasoning follows closely the one of 
Ref. [g8]. 
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Table 5. All Z2 topological insulators (superconductors) can be described through 
dimensional reduction as lower dimensional descendants of parent Z topological 
insulators (superconductor) in the same symmetry class, which are characterized by 
non-trivial winding/Chern numbers. Symmetry classes in which chiral ( "sublattice" ) 
symmetry is present (S=l) are colored blue (the Cartan label marked in bold- face), 
whereas those where it is absent (8=0) are colored red. 



4-1. Topological insulators with chiral symmetry 

Four among the eight symmetry classes of table |] (see also table jT]) are invariant under 
chiral ( "sublattice" ) symmetry, which is a combination of particle-hole and time-reversal 
symmetry. These four symmetry classes are called BDI, CI, CII, and DHL It is possible 
to characterize the topological properties of chiral symmetric Z topological insulators 
by a winding number (topological invariant), which is defined in terms of the block-off 
diagonal projector (see section p. 2.1 ). From this topological invariant we can derive 



a Z2 classification by imposing the constraint of additional discrete symmetries (such 
as C and T) for the lower-dimensional descendants. For a d- dimensional Z topological 
insulator (superconductor) in a given symmetry class, we distinguish between first and 
second descendants, whose (reduced) dimensionality is (i — 1 and d — 2, respectively. 



4- .1.1. Symmetry properties of winding number. In order to better understand the role 
played by T and C, we first study the transformation properties of the winding number 
( ^) and the winding number density (^if ) under these discrete symmetries. First of all, 
we note that the winding number density is purely real, ii'2n-i-ib] — '^'in+iWl^ which can 
be checked by direct calculation. Without any additional discrete symmetries, i.e., for 
symmetry class AIII, the off-diagonal projector g, and hence the winding number density 
(0) are not subject to additional constraints. However, for the symmetry classes BDI, 
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Dili, CII, and CI, the presence of time reversal and particle-hole symmetries relates 
the projector q at wavevector k to the one at wavevector —k. As a consequence, the 
configurations of the winding number density w in momentum space are restricted 
by time-reversal and particle-hole symmetries. Let us now derive these symmetry 
constraints on the winding number density ii'2n.+i[Q'] for any odd spatial dimension 
d = 2n + l. 

Classes Dili and CI: First we consider symmetry classes Dili and CI, where the block 



off-diagonal projector satisfies [26 



q\-k) = eq{kl e={ ^\ ™ (127) 



We introduce the auxiliary functions 
v^{k) : 



(g^^g) (k) and v,{k) := (g^?^) (k)- (128) 



From the symmetry property (|127|) it follows that v^{k) = —v^{—k), irrespective of the 
value taken by e. Noting that (dg"'')g = — gMg we find 



■°-+Hr [v^^{-k)v^^{-k) ■ ■ ■v^.^^^.i-k)]* . (129) 



Thus, the winding number density for symmetry class Dili and CI is subject to the 
constraint 

^2n,+l m] = (-1)"+^^L+1 ki-k)] . (130) 

Consequently, the winding number z/4„4.i[g] is vanishing in d = 4n -|- 1 dimensions, i.e., 
there exists no non-trivial topological state characterized by an integer winding number 
in {4n + l)-dimensional systems belonging to symmetry class Dili or CI (cf. table 
Conversely, in d = {An + 3) dimensions both in class Dili and CI there are topologically 
non-trivial states {ii or iv in table ^), and one would naively expect that for each of 
these states there are lower dimensional descendants characterized by Z2 topological 
invariants. However, this is not the case, as we will explain below. 

Classes BDI and CII: Next we consider symmetry classes BDI and CII, where the 
projector satisfies [^] 

Gq*{-k)G"' = q{k), G=\ l^f/ (131) 

Introducing the auxiliary function v^{k), equation ( |128| ), as before, we find that 
'^fii.k) = —Gv*^{—k)G~^, for both class BDI and CII. From the cyclic property of the 
trace it follows that 

= e-^----+H-l)'"^'tr [v^^[-k)v^,{-k) ■ ..v^,^,,{-k)Y . (132) 
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Figure 4. (a) One-parameter interpolation of two topological insulators represented 
by the projectors qi and q2- In (b), the one-way interpolation (a) is extended by 
making use of the discrete symmetry, (c) Two different interpolations, each colored 
red and blue, respectively. In (d) the original interpolations are rearranged into two 
different interpolations. 



Hence, the winding number density in class BDI and CII is constrained by 

W2n+i [qm = [q{-k)] . (133) 

As a result, the winding number z>'4n~i[q'] is vanishing ind = 4n — 1 dimensions, i.e., there 
are no non-trivial (4n — l)-dimensional Z topological insulators belonging to symmetry 
class BDI or CII (see table On the other hand, in (4n -|- 1) dimensions there exist 
Z topological states in both class BDI and CII {ii or iv in table and one would 
expect, as before, that for each of these states there are lower- dimensional topologically 
non-trivial descendants. Again, this expectation turns out to be incorrect, as we will 
explain below. 

4- 1.2. Z2 classification of first descendants. In this section, we show how a Z2 
topological classification is obtained for 2n-dimensional, chiral symmetric insulators 
(supercondcutors) under the constraint of additional discrete symmetries (i.e., C and 
T). To uncover the Z2 topological characteristics of chiral symmetric systems, we study 
the topological distinctions among the block off-diagonal projectors q{k). Let us consider 
two projectors qi{k) and q2{k), whose momentum space configuration is restricted by 
the symmetry constraint of one of the classes DIII/CI, equation ( |127| ), or BDI/CII, 



equation ( |131| ). We introduce a continuous interpolation q{k, t), t E [0, vr] between these 
two projectors (see figure ^) with 

q(k, t = 0) = qi{k) and q{k, t = n) = q2{k). (134) 

Since the topological space of 2n- dimensional class AIII insulators is simply connected, 
the continuous deformation q{k,t) is well defined. In general, q{k,t) does not satisfy 
the symmetry constraints encoded by equation (|127|) or equation ( |131j ), respectively. 
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However, by combining equation ( p.34| ) with its symmetry transformed partner, we can 
construct an interpolation that obeys the discrete symmetries of the given symmetry 
class (see figure ^p). In the case of symmetry class DIII/CI we define for t G [tt, 27r] 

q{k,t) = eq\-k,2n-t), ^ = | ^J' ™ (135) 

while, in the case of symmetry class BDI/CII we set for t G [tt, 27r] 

q{k,t) = Gq*{-k,2n-t)G-\ G = \ "^'^ J^J^/ (136) 

y lay, Lll. 

Equations ( [134| )-( p!SBD with the parameter t replaced by the wavevector component /c2n+i 
represent a (2n+l)-dimensional projection operator respecting the symmetry constraints 
of the corresponding symmetry class. Consequently, a winding number for q{k,t), 
h'2n+i[q{k,t)], can be defined in the {k,t) space. Two different interpolations q{k,t) 
and q'{k,t) generally give different winding numbers, i^2n+i[(l{k,t)] ^ i^2n+ib'(^5 ^)]- 
However, we can show that symmetry constraint (|135|) or (|136|) , respectively, leads 
to 

z/2„+i [g(fc, t)\ - V2n+i Wik, t)] = mod 2. (137) 



To prove equation ( |137| ) we introduce two new interpolations ri{k,t) and r2{k,t) that 
transform into each other under the discrete symmetry operations (see figure ^) 

q{k,t), tG[0,7r] 
q'{k,27r-t), t G [vr, 27r] , 



ri{k,t) 



These are recombinations of the deformations q{k,t) and q'{k,t) with 

i^2n+i[q{k,t)] - U2n+i[q'{k,t)] = U2n+i[ri{k,t)] + U2n+l[r2{k,t)]. (139) 

Now, we make use of the result from section [4.1. 1| . Namely, for symmetry class BDI/CH 
we found 

w,n+i[q{k,t)] = wl^^M-k, -t)] (140) 

in d = An + 1 spatial dimensions, whereas for symmetry class DIII/CI 

w^n-i [q{k, t)] = wl^_^ [qi-k, -t)] (141) 

in (i = 4n — 1 spatial dimensions. Let us first consider class BDI/CII in d = 4r;, + 1 
dimensions. With equation (20) and ( |140| ) we obtain 

^4n+i[ri]= / d'^"/cdtw4„+Jri(A;,t)] 



J d'-kdtwl^^,[n{-k,-t)] 

J d^''kdtwl^^,[r2{k,t)] = i^4n+i[r2], (142) 
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where we made use of the fact that ri transforms into r2 under the discrete symmetry 
operations of class BDI/ CII. In conclusion, we have shown v^n+i [9(^5 i)\—^in+i [(I'ik, t)] = 
'^J^4n+i[ri{k,t)] G 2Z for any two interpolations q{k,t) and q'{k,t) belonging to class 
BDI/CII. Hence, we can define a relative invariant for the 4n-dimensional projection 
operators qi{k) and q2{k) 

'^An[qi{k),q2{k)] = {-ir-^'^'^^'''^\ (143) 

which is independent of the particular choice of the interpolation q{k,t) between qi{k) 
and q2{k). Once we have identified a "vacuum" projection operator, e.g., qo{k) = go, 
we can construct with equation (|143|) a Z2 invariant: non-trivial Hamiltonians are 
characterized by z/4„[g, go] = —1, whereas trivial ones satisfy z/4„[g, go] = +1. Finally, we 
note that the calculation leading to the relative invariant ^'4„,[gi, g2], equation (|143|) , can 
be repeated for symmetry class DIII/CI in c? = 4n — 1 dimensions yielding a relative 
invariant i'4n-2[qi, q2]- 

For d = 4n — 1 {d = An + 1) the dimensional reduction seems to be possible both 
for class Dili and CI (BDI and CII). However, for a given n the dimensional reduction 
seems to be meaningful only for one of them; while for one of them (corresponding to 
a in table ^), there is a decendent Z2 insulator in one dimension lower, for the other 
(corresponding to iv in table ^), there are no decendent Z2 insulators. In other words, 
for a given n = odd (even), there are lower dimensional Z2 topological insulators for 
either one of class Dili and CI (BDI and CII), but not both. This is because the above 
procedure does not apply when the classification of topological insulator is not Z, but 
2Z. 



4.1.3. Z2 classification of second descendants. The dimensional reduction procedure 
presented in the previous subsection can be repeated once more to obtain a Z2 
classification of the second descendants. As before, we first focus on symmetry class 
BDI/CII and study the topological distinctions among the block off-diagonal projectors. 
We consider two (4n — l)-dimensional projectors gi(/c) and q2{k), which satisfy the 
symmetry constraints imposed by class BDI/CII. We define an adiabatic interpolation 
g(A;,t),t e [0,27r] 

q(k, t = 0) = gi(/c), q{k, t = n) = g2(/c). 




q{k,t) = Gq*{-k,-t)G-\ G={::' ™. (144) 

We can interpret q{k,t) as a 4n-dimensional projector belonging to symmetry class 
BDI/CII. Therefore, for two deformations q{k,t) and q'{k,t) of the form ( |144 ) a 
relative invariant h'in[q{k,t),q'{k,t)] can be defined, as discussed in section |4.1.2| . It 
turns out that due to condition (|144|) the invariant h'^n[q{k,t),q'{k,t)] is independent 
of the particular choice of interpolations, i.e., i/4„[g(/c, t), g'(/c, t)] = +1 for any two 
deformations q{k,t) and q'{k,t) satisfying condition (|144|) . In order to prove this, we 
consider a continuous interpolation r{k,t,s) between the two deformations q{k,t) and 
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q'{k, t) with 

r{k, t, s = 0) = q{k, t), r{k, t,s = tt) = q'{k, t), 
r{k,t = 0,s) = qi{k), r{k,t = vr, s) = q2{k), 

r{k,t,s) = Gr*{-k,-t,-s)G'\ G = | I""' J^J^/ (145) 

1^ lay, bll. 

This represents a (4n+ l)-dimensional off-diagonal projector in symmetry class BDI/CII 
with the winding number i/4„+i [r(/c, t, s)]. We note that r{k,t,s) is not only a 
deformation between q{k,t) and q'{k,t), but can also be viewed as a continuous 
interpolation between r{k,0,s) = qi{k) and r(A;,7r, s) = q2{k) (for any s G [0, 27r]). 
Therefore, we find u^nlqik, t), q'{k, t)] = iy^n[r{k, t, 0),r{k, t, vr)] = u^nirik, 0, s),r{k, vr, s)] 
Since r(/c,0,s) = qi{k) and r(/c, 7r,s) = g2(^) are independent of s, we find 
that z/4„[r(/c, 0, s), r(/c, vr, s)] = (— l)'^4n+iH _ _|_]^^ Hence, we have shown that 
^in[q{k,t),q'{k,t)] only depends on qi{k) and g2(^)- Therefore, i^4n[qo,q{k,t)] together 
with a reference ("vacuum") projector go constitutes a well-defined Z2 invariant in 4n — 1 
dimensions. 



4-1 .4- Example: (i = 3— 7'2— t-I. As an example we consider a 3-dimensional 
topological Dirac superconductor belonging to symmetry class DHL Consider a. d = 3 
dimensional Dirac Hamiltonian, 

'H^gy^(fc, m) = k^a^ + kyUy + /c^a^ — mi/37^. (146) 

This is nothing but the chiral topological Dirac superconductor in class DHL This 
Hamiltonian is essentially identical to the BdG Hamiltonian describing the Bogoliubov 
quasiparticles in the B phase of superfluid ^He ||2^, |70|, ^ ||. It also describes an auxiliary 



Majorana hopping problem for an interacting bosonic model on the diamond lattice . 
In this basis, discrete symmetries are given by 

r : (ay^ r,) [^fgf (-/c, m)] * (a^, ® r,) = Up^ik^m), 

C : {ay ® Ty) [n1^\-k,m)X {ay ® Ty) = -nfff{k,m). (147) 

Combining these two, we have chiral symmetry, which allows us to define the winding 



number 1/3. From the calculations in Sec. |2.5| , the winding number is non-zero, 
= ±1/2, depending on the sign of the mass. 



By dimensional reduction, we obtain the Hamiltonian in one-dimension lower 

^?5) 



'Hf^)'^{k, m) = k^a^ + kyUy - mi(3'j^. (148) 

This Hamiltonian is a 2D analogue of "^He-B, and is unitarily equivalent to the direct 
product of spinless p + ip and p — ip wave superconductors. As this is obtained from 

jj In Rcf. [55( 1, a topological invariant counting the number of gap-closing Dirac points in an extended, 
higher-dimensional parameter space was discussed for a four-band example ( |146| ) . as opposed to 1/3 [q\ 
defined for a given topological phase. Moreover, the crucial role for the protection of topological 
properties arising from the combined time-reversal and charge- conjugation symmetries of symmetry 
class Dili was not discussed. 



CONTENTS 



47 



dimensional reduction of a parent topological superconductor ( p.46| ) this is a Z2 state. 
Finally, by further reducing dimensions, 



kxUx - mi/37' 



(149) 



This is a ID ^^.-wave superconductor. Again, this is a Z2 state. 

We now discuss the topological character of these states in more detail. The Z2 
nature of the state in this d = 2 example can be studied by the Kane-Mele invariant 



MM 



It can be expressed as an SU(2) Wilson loop 
1 



Wsv{2)[L] := -trPexp 



Aik) 



(150) 



Here A{k) is the SU(2) Berry connection, P represents the path-ordering. By definition, 
W^su(2)[-Z^] is a well-defined and gauge invariant quantity for any loop L in the BZ. For 
time-reversal invariant systems, it is useful to consider a loop L that satisfies time- 
reversal symmetry; i.e., a loop which is mapped onto itself (up to reparameterization) 
under k — )■ —k. For these time-reversal invariant loops, the SU(2) Wilson loop is 
quantized. 



Wsvi2)[L]= ±1, 



:i5ii 



and provides a way to distinguish different Z2 states [T5|, ^ . The quantization of the 
Wilson loop can be proved by first noting that because of T, Bloch wavefunctions at k 
and —k are related to each other by a unitary transformation w{k), 

\Qulik))=w,f,ik)\u7{-k)) (152) 
where represents T operation and 

^dk) ■■= {u7i-k)\Qu7{k)). (153) 
Accordingly, the Berry connection at k is gauge equivalent to Aj^ at —k, 

A^{-k) = + w{k)Al{k)w\k) - w{k)df,w\k) 

= -w{k)A*^{k)w\k) -w{k)df,w\k). (154) 

Because of this sewing condition ( |154|) of the gauge field, when plugged in (|150|) , 
contributions at k and —k in the path integral cancel pairwise, except at those momenta 
which are invariant under T by themselves. Thus, 

Wsvi2)[L] = l[Pf[w{K)], (155) 

K 

where K is a. momentum which is invariant under k — )■ —k. 

We now compute the Z2 number for the d = 2 class Dill topological superconductor 
( |148|) . From (|8|), 

/ —kx + iky \ / im \ 



\ulik)) 



V2\ 



V 



im 

A 



1^2 (fc)) 



V2X 



im 

'kx ^ky 
X 




(156) 
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Then, their time-reversal counterparts are 

/ A \ 



\Qu~,{k)) 



V2X 





— im 

\y kx iky J 



\eu^{k)) 



V2X 



( 



\ 





-A 

~kx ~\~ iky 
im 



(157) 



where G = [cTy ® Tx)}C with /C being the complex conjugation. The "sewing matrix" 
w{k) can then be computed as 



w{k) := {Uai-k)\euf^{k)) 



1 



ikx + kjj 



m{k) 



(15^ 



A(A;) I —m{k) +ikx + ky 

Here, we have regularized Dirac Hamiltonian properly, by making the mass /c-dependent, 
e.g.,. m — )■ m{k) = m — Ck"^. With this regularization, k = and k = oo are the two 
time-reversal invariant momenta. One finds 

m(0). 



m 



rPf (ia 



Pf^(oo) = ?^Pf (ia„ 



A(oo) 



(159) 



and hence 

Wsvi2)[L] = -sign(m)sign(C). (160) 

That is, when sign(m) = sign(C) the Hamiltonian (|146|) is a Z2 topological 
superconductor. 



Topological insulators lacking chiral symmetry 

Four among the eight ( "real" ) symmetry classes of table |^ break chiral ( "sublattice" ) 
symmetry. The Z topological insulators (superconductors) which break chiral symmetry 
are characterized by a Chern number [see equation ([Tl|) ]. Using this topological 
invariant one can derive the Z2 classification for the lower dimensional descendants. This 



derivation is analogous to the discussion in section 4.1 and has been performed previously 



in Ref. p8[ for the Z2 topological insulators which break chiral symmetry. For these 



reasons we do not repeat the argument here and refer the reader to Ref. p8[ for details. 
For the Z2 topological insulators (superconductors) which break chiral symmetry, one 
can construct a Z,2 index similar to the one constructed by Moore and Balents [O, in 



all (even) dimensions 74 



5. Discussion 

In this paper we have performed an exhaustive study of all topological insulators and 
superconductors in arbitrary dimensions. The main part of this paper deals with 
dimensional reduction procedures which relate topological insulators (superconductors) 
in different dimensions and symmetry classes. We also discussed topological field 
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theories in D = d+ 1 space time dimensions describing linear responses of topological 
insulators (superconductors). Furthermore, we studied how the presence of inversion 
symmetry modifies the classification of topological insulators (superconductors) (see 
appendix |C|). In the following we give a brief summary of the main results of the paper. 



5.1. Dimensional reduction procedures 

We have constructed for all five symmetry classes of topological insulators or 
superconductors a Dirac Hamiltonian representative, in all spatial dimensions. Using 
these Dirac Hamiltonians as canonical examples, we have demonstrated that topological 
insulators (superconductors) in different spatial dimensions and symmetry classes can be 
related to each other by dimensional reduction procedures. These Dirac representatives 
have been useful for constructing string theory realization of topological insulators 
and superconductors: In Ref. [7S], a one-to-one correspondence between the ten-fold 



classification of topological insulators and superconductors and a K-theory classification 
of D-branes was established, where open string excitations between two D-branes of 
various dimensions are shown to reproduce all Dirac representatives constructed in this 
paper. 



5.1.1. The "complex" case. We first studied topological insulators (superconductors) 
with Hamiltonians which are complex (i.e., belonging to the "complex case" of table 
that is, classes A and AIII from table |l]). Here, starting from a topological insulator 
which lacks chiral symmetry (i.e., class A) in even spatial dimensions d = 2n, we 
obtain a topological insulator (superconductor) in d = 2n — 1 which possesses chiral 
symmetry (i.e., class AIII), by dimensional reduction. In order to understand how 
the topological characteristics of the d = 2n — 1 class AIII topological insulator 
(characterized by the winding number t'2n-i) are inherited from its higher dimensional 
parent, the d = 2n class A topological insulator (characterized by the Chern number 
Ch„), it is important to first point out a few properties of the momentum space topology 
of Bloch wavefunctions in symmetry classes A and AIII. First we note that the non-zero 
Chern number in class A leads to an obstruction to the existence of globally defined 
Bloch eigenfunctions. I.e., it is not possible to construct Bloch eigenf unctions for the 
parent class A topological insulator globally on the d = 2n dimensional Brillouin zone. 
Hence, the Bloch eigenfunctions can only be defined locally on some suitably chosen 
coordinate patches. For the descendant topological insulator (superconductor) with 
chiral symmetry, on the other hand, since there is no Chern invariant ind = 2n — 1, there 
always exists a basis (i.e., a gauge) in which the Bloch eigenfunctions are well-defined 
globally on the entire d = 2n — 1 dimensional Brillouin zone. Such a basis is provided 
by the one, in which the Hamiltonian is block off-diagonal. The existence of such a 
block off-diagonal basis, in turn, is implied by the chiral symmetry of the descendant 
class AIII topological insulator (superconductor). As a consequence, the role played by 
the chiral symmetry is to guarantee that Bloch eigenfunctions of class AIII topological 
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insulators (superconductors) can be defined globally on the entire Brillouin zone. 

Now, by adopting such a block off-diagonal basis for the class AIII topological 
insulator, we have shown that the d = 2n — 1 dimensional Brillouin zone of the 
descendant class AIII Hamiltonian defines the boundary of two local coordinate patches 
of the d = 2n dimensional Brillouin zone of the parent class A topological insulator. 
The transition function between these two coordinate patches is given by the block off- 
diagonal projectoiU] q of the descendant Hamiltonian, and the Chern number Ch„ of 
the parent Hamiltonian can be written in terms of a winding number of this transition 
function (see section |27^) . 



5.1.2. The "real" case. Also demonstrated is, using the Dirac Hamiltonian 
representatives, the periodicity 8 shift of symmetry classes for the topological insulators 
(superconductors) in the 8 symmetry classes in which the Hamiltonian possesses at least 
one reality condition (arising from T or C). Once one has constructed a representative 
in terms of a Dirac Hamiltonian, the reality properties of the spinor representations 
of the orthogonal groups SO(A^) (which are linked to the reality properties of the 
representations of the Clifford algebra formed by the gamma matrices) leads directly 
to this 8- fold periodicty in the spatial dimension d (see section From the Dirac 
Hamiltonian representatives of topological insulators (superconductors) in all d spatial 
dimensions (all of which are of course massive), one can realize a. d — 1-dimensional 
boundary, by making a domain wall at which the mass term changes sign [j7U|. There 
then appears a boundary state localized at the boundary, which is a defining property 
of topological character in the bulk. These boundary states are of special kind, as 
they are protected from the appearance of a gap and from Anderson localization, by 
the time-reversal and charge- conjugation symmetry properties of the specific symmetry 
class. Indeed, one of the underlying strategies of classification adopted in Ref. 



was to classify the properties of these boundary Hamiltonians in each symmetry class, 
including those of Dirac type, which must be completely gapless and delocalized when 
they come from a topological bulk, but can otherwise be gapped and localized for a 
non-topological bulk. By constructing massive Dirac Hamiltonian representatives in 
the bulk for all topological insulators (superconductors), which give rise to the gapless 
boundary Dirac Hamiltonians, we undertook here in some sense a complementary task 



as compared to that in Ref. [26 



5.2. Topological field theories 

In the low-energy limit, all topological insulators (superconductors) can be described 
by topological field theories for the linear responses in space-time, which characterize 
universal (in principle) experimentally accessible observables of the topological features, 
such as, e.g., in transport. For example, for symmetry class A in even spatial dimensions 
d = 2n, the topological field theory of the linear responses in D = d + 1 space-time 

t of the type described in equation (34) of the first article in Ref. 
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dimensions is given by the Chern-Simons action, whose couphng-coefficient is the nth 
Chern number Ch„. For symmetry class AIII in odd spatial dimensions ci = 2n — 1, on 
the other hand, the topological field theory m D = 2n space-time dimensions is given 
by the 9 term, whose coefficient (the 9 angle) is given by the winding number f2n+i 
(see section p^ ). For topological singlet superconductors (symmetry classes C and 
CI), one can couple external SU(2) gauge field to the conserved spin current operator 
of the BdG quasiparticles. The spin response in topological singlet superconductors, 
can be described by SU(2) gauge theory with Chern-Simons type topological term in 
c? = 2n [§] and with the ^-type-term in c? = 2?t, + 1 |7^. In passing we note, that it is 
also possible to establish a connection among these topological field theories in various 
space-time dimensions and symmetry classes via dimensional reduction procedures. The 
topological field theory formulation may be a good starting point to explore, more 
generally, topological phases in interacting systems beyond those that are presently 
known. 



5.3. Topological insulators with inversion symmetry and either time reversal or charge 
conjugation symmetry 

We have also studied the restrictions imposed on groundstate properties of topological 
insulators (superconductors) by the presence of inversion symmetry. That is, we studied 
topological states that are protected by a combination of spatial inversion (denoted by 
X) and an additional discrete symmetry (i.e., T or C). These systems are invariant under 
the combined symmetry operations T ■ X or C ■ X, but all three symmetries T, C and X, 
are assumed to be absent. We have determined the space of projectors describing these 
topological states. From the homotopy groups of these space of projectors follows the 
classification of topological states that are protected by T ■ X or C ■ X (see table |C^ ). 



5.4. Directions for future work 

An important direction for future study is the search for experimental realizations 
of three-dimensional topological singlet or triplet superconductors. Given how fast 
experimental realizations of the QSHE m d = 2 and the Z2 topological insulators 
in (i = 3 have been found, we anticipate a similar development for the three- 
dimensional topological singlet or triplet superconductors. For example, unconventional 
superconductors in heavy fermion systems, typically possessing strong spin-orbit effects, 
have been studied extensively over the years. They might be good candidates for 
topological superconductors. Moreover, the search for non-trivial topological quantum 
ground states is not restricted to free fermions, or BCS quasiparticles, but includes 
also, more generally, strongly interacting systems other than BCS with emergent free 
fermion behavior at low energies. Furthermore, it should be noted that the classification 
scheme given by table |]is also applicable outside the realm of condensed matter physics. 
For example, topological properties of color superconducting phases [^, |7^, which are 



predicted to occur in quark matter, can be discussed in terms of the present classification 
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scheme. 
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A. Cartan symmetric spaces: generic Hamiltonians, NLcrM field theories, 
and classifying spaces of i^-theory 

In this appendix we review the appearance of Cartan's ten-fold list of symmetric spaces 
in the context of (i) basic quantum mechanics: where they describe the time evolution 
operators exp(it'H) of generic Hamiltonians 7/, (ii) NLcrM field theories: where they 
describe the "target space" manifold of the NLcxM, and (iii) K-Theory: where they 



describe the "classifying space" (briefly discussed at the end of subsection of the 
introduction). As reviewed in the introduction, the homotopy groups of these symmetric 
spaces play a key role in the classification of topological insulators and superconductors 
both in the approach of Ref. [EH] , which makes use of results from Anderson localization 



physics, as well as in the approach of Ref. [B^, which is based on K-theory. 



Cartan 


time evolution op. 


fermionic replica 


classifying 


label 


exp{it?^} 


NLcrM target space 


space 


A 


U(iV) X \]{N)/\]{N) 


U(2n)/U(n) x U(n) 


U(iV + M)/U(iV) X U(Af) = Co 


AIII 


\]{N + M)/\]{N) X U(M) 


U(n) X U(n)/U(n) 


U(iV) X U(iV)/U(A^) = Ci 


AI 


U(A^)/0(iV) 


Sp(2n)/Sp(n) x Sp(n) 


0(iV + M)/0{N) X 0(M) = Ro 


BDI 


0{N + M)/0{N) X 0(M) 


U(2n)/Sp(2n) 


0{N) X 0{N)/0{N) = Ri 


D 


0(iV) X 0(iV)/0(A^) 


0(2n)/U(n) 


0(2iV)/U(iV) = i?2 


Dili 


S0(2iV)/U(iV) 


0(n) X 0(n)/0(n) 


V{2N)/Sp{2N) = i?3 


All 


U(2iV)/Sp(2iV) 


0(2n)/0(n) x 0(n) 


Sp{N + M)/Sp(A^) X Sp(M) = R4 


CII 


Sp(A^ + M)/Sp(A^) X Sp(M) 


U(n)/0(n) 


Sp(iV) X Sp{N)/Sp{N) = R5 


c 


Sp(2iV) X Sp(2iV)/Sp(2iV) 


Sp(2n)/U(n) 


Sp(2iV)/U(A^) = Re 


CI 


Sp(2Af)/U(Af) 


Sp(2n) X Sp(2n)/Sp(2n) 


U(iV)/0(iV) = Rr 



Table Al. Three appearances of the hst of Cartan's ten symmetric spaces. First 
column: unitary time evolution operator, second column: (compact) target space of 
NLcrMs, third column: classifying space. (We use the convention in which m = even 
in Sp(to). - Moreover, the Cartan labels of those symmetry classes invariant under the 
chiral symmetry operation S — T -C from table || are indicated by bold face letters.) 



In table ^T] we list for each symmetry class denoted by its Cartan label in the 
first column, the time evolution operator in this symmetry class (penultimate column 
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of table |T]) in the second column, the target spaces ||| of the NLaM field theories in the 
third column, and in the last column the classifying space appearing in K-Theory [ p7[] . 

Here we would like to re-emphasize the remarkable fact that the same ten Cartan 
symmetric spaces describe all three objects listed in the three last column of table [M . 
Furthermore, there are remarkable relations between these columns. First, the second 
column ("time evolution op.") can be obtained from the last column ("classifying 
space") by shifting the entries in the last column down by one entry (modulo eight, and 
modulo two in the "real" and "complex" cases, respectively). Second, the third column 
("NLcrM target space") is obtained from the last column by performing a reflection 
(modulo 8) in the last column about the entry in the row labeled D (and an exchange 
for the complex case). Consequently, there is then also a resulting relationship between 
the second and the third columns . 

As reviewed in subsection |1.2| of the introduction, the classifying spaces listed in 
the last column of table ^1] describe topological insulators (superconductors) in zero 
dimensions, i.e., at one point in space|§[ The disconnected components of this space of 
Hamiltonians Q (which cannot be continuously deformed into each other) are labeled 
by Z or Z2 appearing in the = column of table ^, which denotes the list of zeroth 
homotopy groups of the spaces in the last column of table |A1| . All higher homotopy 
groups of the same spaces can then be inferred from table |^ due to the dimensional 
periodicity and shift properties visible in that table. 



B. Spinor representations of SO(A^) 



In this appendix we review spinor representations of SO(A^) and their properties under 
reality conditions [81, 32|. It is most convenient to discuss spinor representations in 
terms of Clifford algebras, i.e., in terms of gamma matrices {r^^^}a=i,,,,,Ar satisfying 
{r|^^, r^^^} = 25jk, with j,k = 1, . . . ,N. Given such a set of gamma matrices a spinor 
representation of SO(A^) can be readily obtained 

i 



M. 



jk 



(B.l) 



with the SO(A^) generators M. 



jk- 



B.l. Spinors of SO {2n + 1) 

In what follows we will focus on SO(2?7, + 1), which has a 2 "-dimensional irreducible 
spinor representation. The gamma matrices in the Dirac representation for iV = 2?7, + 1 

I There are three different varieties of NLcrM field theories: supersymmetric, fermionic replica, and 
bosonic replica models. In the latter two formulations, the replica limit must be taken, where one lets the 
number of replicas, N, tend to zero at the end of the calculations. While the target spaces of fermionic 
replica models are compact, those of bosonic replica models arc non-compact. In the supersymmetric 
formulation, on the other hand, the target spaces are supermanifolds with both bosonic and fermionic 
coordinates. The fermionic and bosonic replica NLcrMs, with N finite, can be viewed as fermion-fermion 
and boson-boson subsectors, respectively, of the corresponding supersymmetric NLg Ms Ip8| . 
§ Ref. compare also table III (second column, with fc ~ 0) of the first article of |26|]. 
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are defined recursively by 



r(2n+i) = r^2n-i) ® ^3, a = 1, ■ ■ ■ , 2n - 2, 



p2n+l / ;\npl p2 p2n m 

^ f2n+l) - ^ (2n+l)J- (2n+l) (2n+l)' l^'^J 



(2n+l) 

^2n+l _ 
(2n+l) — \ '■I {2n+l)J- (2n+l) ' ' ' {2n+l)' 

where Jon-i is the 2""^ X 2"~i identity matrix. [The Eramma matrices in the Dirac 
representation for N = 2n can be constructed by just leaving out ^'^(^^+1)^ 
^{2n) ~ ^'{2n+i)^ with a = 1, ■ ■ ■ , 2n.] To be more explicit, 

r(2n+l) = CTl ® ^3 (8) ■ ■ ■ (g) (73 , 
n-1 

r(2n+l) = Cr2 ® 0-3 O ■ ■ • (g) 0-3, 

^3 
(2n+l) 

^4 
(2n+l) 



n-1 



rf2n+ll = O-Q ® ® 0-3 ® ■ ■ ■ ® CTS, 



n-2 



rf2n+l1 = CTQ ® 0-2 ® 0-3 (g) ■ ■ ■ (g) 0-3, 



n-2 



(2n+l) 
^2n 



n-1 



r'aVi) = 5^o®- - -^c7o ®(T2, (B.3) 



n-1 



and 



r?2n+l) = .^3®---®a3. (B.4) 



From the explicit construction of the gamma matrices, we infer that r(2n+i) 
all real, and r(2n+i)" purely imaginary. In order to implement discrete symmetries 
on the space of Dirac Hamiltonians we define the matrices 

rjl pi p3 p2n— 1 

-°(2n+l) •— (2n+l)J- (2n+l) ' " >- (2n+l)' 

-^(^2n+l) r(2„+i)r^2n+l) ' ' ' r^2n+l)- (B-5) 

By use of the reality and anti-commutation properties of the gamma matrices, one finds 
that 

[5(Wl)]*42n+l) = (-lr^"-^^/^ 

[5?2n+l)]*i^(Wl) =(-1)"^"^^^/^- (B.6) 

The operator -B(^2n+i) used to construct Majorana (real) representations of S0(2n-|-1). 
For odd = 2n + 1, they are possible when ||8l], ^ 

[5(Wi)]*i?(Wi) = (-ir^"^^^/^ = i 

^ n = 0,3 mod 4 {N = 2n + 1 = 1,7 mod 8) . (B.7) 
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Observe that when n = even all gamma matrices can be made real and symmetric, 
whereas when n = odd, they can be made purely imaginary and skew-symmetric. 
For later use, we also introduce 



r12 rR2 rI — ( i\n-p 

-^{2n+l) [^{2n+l)\ ^{2n+l) " ^ 



2n+l 
(2n+l)' 



r1 p2n 
^(2n+l)'- (2n+l)5 



'{2n+l) ■ 

rS2 ry2 -p2n 

-°(2n+l) -"(2n+l)-'- (2n+l)' 

Note that these matrices satisfy 



(B.8) 



Bl2n+l)^'(2n+l) [^{2n+l)] 



-^(271+1)^(271+1) [-^("2n+l)] 



r12 -pa [r>12 1 

^(2n+l)^ (2n+l) L-D(2n+l)J 



(' 1 Nn+lpa* 

I ^) J-(2n+l)' 

y ^) (2n+l)' 
"J-J J- (2„+l), " - 

~-'^"2n+l)' ^ = 
+^"271+1)' 



a = 1, ■ 

a = 2n 



,2n, 
1, 



= 1,- ■■,2n + 1, 

l,---,2n, 
2n + l. 



(B.9) 



In the following subsection we will use the "5- matrices" , equations (p.5| ) and ( p.8| ), 
as symmetry operators in order to implement discrete symmetries on the space of Dirac 
Hamiltonians. To identify the character of these discrete symmetries we first compute 
the sign r]B^2n+i) picked up by 5(2„+i) under transposition, [5(2„+i)]'^ = '7B(2„+i)5(27i+i)- 
Note that the sign ?7S(2„+i) independent on the choice of basis for the gamma 
matrices (see Ref. From equations ( [B.6| ) and ( |B.7| ) and by use of the property 

1, it follows that 

{2n+l) 

Similarly, 



[-^(^2n+l)]^-^(^2n+l) 



(B.IO) 



''°{2n+l) 
'°{2n+l) 
{2n+l) 



(_l)"("-l)/2^ 
_ (_l)"("+l)/2 
. ^x„(„+3)/2^ 



(B.ll) 



B.2. Discrete symmetries of Dirac Hamiltonians 

Let us now determine the symmetry properties of the Dirac Hamiltonians ^^2n+^3)^(^)' 

'^{Tn+3)^(k)^ • • '^?2n?3)^'^^'^)' ^^^^^^ equatlous (|100|) through (^). We find that 
these Dirac Hamiltonians satisfy the following symmetry conditions 



(0) : Bf2n+i)'H'l^2n+l)'i^) [^(Wl)] 



\7l+l 



n 



d=2n+lj 



(B.12) 



'-(27i+iy( ^) 

(0 : Bl,^^,^'Hp::,^{k,m) [i?(Wi)]"' = (-1)" WTn%i-k,m)Y, (B.13) 

(22) : 5f2n+l)^fn?iT(^'"^) [^"+1)] = " 



^?2:+V(-^'"^) 



(2n+l)^(27i+l) 



■ ci=2n-l 



{k, m) 



Bl 



(2n+l) 



= (-1 



\n+l 



5 

1 * 



(B.14) 
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(B.15) 



[IV] 



= {-l 



n+l 



-Sf2n+l)'^(2n+l)^(^' 



^(2n+l) 



(B.16) 



A few remarks are in order 

/d=2n-l /■ u ™^ „„J i/d=2n-3 



The Hamiltonians 'H^^^^^ {k, m) and 'H(^2n+i) satisfy two different discrete 

symmetry conditions, and are therefore also left invariant under the combination 
of these two symmetries, which defines a chiral symmetry. In other words, both 
"^(271+^1) ^ ( ^' ''^) ^^"^ '^t^n+Vj^i^^^) anticommute with a unitary matrix 



^?2:+V(^,^),r?2n-i)| = o. (B.17) 

• The Hamiltonians 'Hf^n+3)ik^^) and ^^2n?3)^(^' 

m) can be made block diagonal 

because 

['^(2n+3)(^5 '"^)' r"2n+3)r(2n+3)] = 0' 

^fn;3T(^' "^)' r(2n,+3)r?2n,+3)] = 0, (B.18) 

where (a, 6) = (2n + 1, 2n + 2), (2r2 + 2, 2n + 3), or {2n + 3, 2n + 1). 

• In the case of the massive Hamiltonian H'f^^^^^^ik^m), the parity transformation!]] 
can be implemented by B^^n+i) 

BlL+i)npn%{k,m) [Sf^l+i)]"' = npn%)i'k,m). (B.19) 

By combining the discrete symmetry -B(2n+i) and the parity symmetry B^^n+i) ' 
Hamiltonian also satisfies 

i^(Wl)^(2=n+l)(^,H W2n,+1)]"' = ("l)'^ [^(2f+"l) ^ • (B.20) 

This discrete symmetry is unique in that, unlike T or C which relates Bloch 
Hamiltonians at k and —k, it constrains the form of the Hamiltonian at a given k. 
This is nothing but the real/pseudo-real condition for S0(2?t,+1). As a consequence, 
the Hamiltonian in fc-space can be written as a real/pseudo-real matrix. Hence, 
classifying topological insulators (superconductors) satisfying the combination of 
B^2n+i) and amounts to classifying real bundles. (It is not necessary 

to satisfy both, though.) An example of a topological insulator /superconductor 
satisfying the combination of B^^n+i) and B^^n+i) constructed in appendix 0. 

II Parity here simply means k — > — fc symmetry. In odd space time dimensions, it is actually not called 
parity. 
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C. Topological insulators protected by a combination of spatial inversion 
and either time-reversal or charge-conjugation symmetry 

In the main text, we have focused on the role of generic symmetries such as time reversal 
and charge conjugation, which are not related to any spatial symmetries. However, 
real systems often have discrete spatial symmetries, such as parity, reflection, discrete 
rotations, etc. Hence, it is meaningful to study the restrictions imposed by these 
symmetries on the Bloch wavefunctions, and thus on the ground state properties. In 
this appendix, we briefly discuss topological states that are protected by a combination 
of spatial inversion and an additional discrete symmetry, such as T or C. 

C.l. Inversion symmetry combined with another discrete symmetry 
Consider a tight-binding Hamiltonian, 

H =^^l,\r)H{r,r')tl){r'), H\r ,r) = H{r,r'), (C.l) 

r,r' 

where if){r) is an n/-component fermion operator, and index r labels the lattice sites. 
(The internal indices are suppressed.) Each block in the single-particle Hamiltonian 
'H[r,r') is a. Uf X Uf matrix, and we assume the total size of the single particle 
Hamiltonian is njV xufV, where V is the total number of lattice sites. The components 
in ipi^r) can describe, e.g., orbitals or spin degrees of freedom, as well as different sites 
within a crystal unit cell centered at r. 

Provided the system has translational symmetry, 

H{r,r') = H{r-r'), (C.2) 

with periodic boundary conditions in each spatial direction (i.e., the system is defined 
on a torus T*^), we can perform the Fourier transformation and obtain in momentum 
space 

H= Y,i^\k)H{k)^{k), (C.3) 

fcGBZ 

where 

^(r) = yy"' ^ e'^^XA;), H(A;)=^ e-*'="H(r). (C.4) 

fcGBZ r 

[Here, if different components in ip[r) were to represent different sites within a unit 
cell centered at r, we could include the phase e^'''^^~^°-^ in the definition of ^{k), 
ipik) — )■ diag (e'^''-''~''°-*)^_-^ ^ '4'{k) where represents a location within a unit cell.] 

In the main text of this paper, we considered topological insulators (superconduc- 
tors) protected by discrete symmetries, such as T and C. We now focus on the effect of 
inversion symmetry. By definition, an inversion Q relates fermion operators at r and at 
-r, 

ni/j{r)n-^ = U'^{-r), (C.5) 
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Table CI. The ten generic symmetry classes of single-particle Hamiltonians classified 
in terms of the presence (Q) or absence (x) of the time-reversal symmetry (ey = +1) 
with T= —1 (rjv = —1) and T= -1-1 {tjv = and the particle- hole symmetry 

(ey = —1) with C= —1 {iiv = —1) and C= +1 {rjv — +1), as well as chiral (or 
sublattice) symmetry denoted by S; see table |^. For the notation of these symmetry 
classes in terms of (ey, rjv), see the main text and Ref. |2^ . 

where [/ is a constant Uf x Uf unitary matrix. (It is important to distinguish between 
inversion and parity symmetry. The former means r — )■ — r for any space-time dimension. 
Parity transformation in even space-time dimension agrees with inversion, but in odd 
space-time dimension, it is different from inversion.) One can imagine, for example, U 
represents a parity eigenvalue of each orbital at site r {s,p, . . . , orbitals), which can be 
either -|-1 or —1. Similarly if different components in iplr) represent different sites in 
the unit cell centered at r, inversion transformation sends r to — r and at the same time 
can reshuffle locations of electron operators within a unit cell, which can be described 
by the unitary matrix U. The invariance of H under Q implies 

U~^n{r,r')U = H{-r,-r'), U-^H{r)U = H{-r) . (C.6) 

Thus, the Bloch Hamiltonians at k and —k are unitarily equivalent, 

U~^V,{k)U = ni-k). (C.7) 

To summarize, definition ( |C.5| ) implies that (i) the spectrum at k is identical to the one 
at —k, and (ii) the unitary matrix U relating the Hamiltonian at k to the one at — k 
is /c-independent. This form of symmetry does not necessarily correspond to inversion 
symmetry verbatim, but might describe an inversion symmetry which is realized as a 
projective symmetry (e.g., invariance of the Hamiltonian under inversion followed by 
a gauge transformation). The only big assumption here is that we have assumed U is 
/c-independent. 

We now combine the inversion symmetry with another discrete symmetry, 
V-^H{k)V = evUi-kf , ey = ±1, 
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with = rjvV, rjy = ±1. (C.8) 

As before we distinguish four different cases [see section and in particular equations 
(§) and (H)]: (ey,?7y) = (1,1) corresponds to a time-reversal symmetry with T= +1, 
{€v,riv) = (1,-1) denotes a time-reversal symmetry with T= —1, {ev,riv) = (—1,1) 
is a particle-hole symmetry with C= +1, and (ey,?7y) = (—1,-1) is a particle-hole 
symmetry with C= —1. These four cases together with the corresponding "Cartan 
label" are summarized in table |Cll . Then, for the combined transformation, 

W-^H{k)W = ewHkf, W := UV, 

with ew = ev- (C.9) 
Iterating this transformation twice, 

H{k) wnikfw-^ w{w-Y'H{k)w^w-\ (C.IO) 

which suggests, from Schur's lemma, W^W~^ oc Hence, as before, for W, we 

distinguish the following two cases, 

W^ = r]wW, w = ±l- (C.ll) 

The signature riw is invariant under unitary transformations. It depends on the 
signature r]v, the signature rju for the inversion U'^ = rjuU, and the commutation 
relation of U and V. Here, note that the signature rju for inversion alone is not invariant 
under unitary transformation, while the signature t]w is invariant. As an illustration, let 
us consider spinless fermions. Time-reversal symmetry T for spinless fermions implies 

H{i,3y =H{i,3)- (C.12) 

Thus, when inversion and T are combined, 

n*{r) = Wn{-r)W-\ H*{k) = W~^n{k)W, where W = U. (C.13) 

In this example and in this basis, the signature of W solely depends on the signature of 
U, rjw = Vu- 

Interesting examples are provided by Dirac Hamiltonians. For a massive Dirac 
Hamiltonian, the mass matrix itself can be taken as the unitary matrix U which was 
introduced in equation (|C.5|) . Put differently, any massive Dirac Hamiltonian possesses 
a U-tjpe symmetry. These U-tjpe symmetries can be, however, different from the 
inversion symmetry which is imposed at the microscopic level. However, if one starts 
from a microscopic lattice model with inversion symmetry, and arrives at a Dirac 
Hamiltonian with a mass in the continuum, inversion symmetry in the continuum must 
be implemented by a mass matrix. 

First, consider Hamiltonian ( p.01| ) 

d=2n 

(^) : ^i2n%ik,m) = J2 ^ar^2n+l) +^rjj;\). (C.14) 

a=l 

This Hamiltonian has a V^-type symmetry with 

^ = i?(Wi), ^v = {-ir, w = r/5i^^^^^ = (-l)"("-i)/^ (C.15) 
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There is inversion symmetry represented by 

U = ^l^:i,y Vu = +l. (C.16) 
The combined transformation W is given by 

V^^ = 5(Wi), ew = {-lT, w = = (-1)"^"^'^/'. (C.17) 

Second, we consider the massive Dirac Hamiltonian ( |103| ) 



d=2n-2 

(m): Hf„2n^-^2(/c,m)= ^ /c.r^^n+i) + '^-M, 

a=l 

With A^:=ir?,^+H_\)r?,Vi)r?,Vi)- (cis) 

The V^-type discrete symmetry for this is 

V = = Vv = VBf,^^,^ = (C.19) 

The inversion symmetry is represented by 

U = M, vu = +1- (C.20) 
Thus, the combined transformation is 

l/K - y\/l±J(2n.+l) - IJ- (2n+l)^ (2n+l)J- (2n+l) ^ (2n+l)J- (2n+l) (2n+l)- 

with = (-1)"+', Vw = (C.21) 
where we have noted UV = VU. 

C.2. Classification of gapped Hamiltonians in the presence of a combination of spatial 
inversion and either time-reversal or charge- conjugation symmetry 

We now classify the ground states of gapped Hamiltonians protected by a W^-type 
symmetry. We distinguish four different cases, {ew^Vw) = (iljil)- 

Let us first consider the case of {ew,Vw) = (1, 1), in which case W is a. symmetric 
unitary matrix, W = W'^ {rjw = +1). Such a matrix is an element of U{nf)/0{nf). Any 
element in U(?7,j-)/0(n/) can be written as W = XX'^ where X is a unitary matrix. 
One can then take a basis in which the Bloch Hamiltonian for any given is a real 
symmetric matrix. In other words, if we define := X^TiX, then "H* = "H. Thus, when 
{ew, Vw) = (+1) +1)) the Hamiltonian is real symmetric and hence Bloch wavefunctions 
can be taken real at each k. We assume there are (A^+) occupied (unoccupied) Bloch 
wavefunctions for each k with A^_|_ + = A'tot(= nf). The spectral projector onto the 
filled Bloch states or the "Q-matrix" in this case can be viewed as an element of the real 
Grassmannian Gn_,n++N- {^) = 0(A^+ + A^_)/[0(A^+) x 0(A^„)]. For a given system, the 
"Q-matrix" defines a map from the BZ onto the real Grassmannian. Hence, classifying 
topological classes of band insulators amounts to counting the number of topologically 
inequivalent mappings from the BZ to Gn_,n++n.(M)- Mathematically, this is given by 
the homotopy group of the space of projectors (i.e., of the real Grassmannian in the 
present case), which can be read off from table 0. 
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1 
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7 ■■■ 


(+1, 


-1) 


Sp(iV + M)/Sp(Ar) X Sp(M) 


z 











Z 


Z2 


Z2 


■■■ 


(-1, 


+1) 


0(2A^)/U(A^) 


Z2 





z 
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Z2 ■■■ 


(+1, 


+1) 


0(A^ + M)/0(A^) X 0(M) 
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Z2 
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■■■ 


(-1, 


-1) 
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Z2 


Z2 
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■ ■■ 



Table C2. Classification of topological insulators (superconductors) protected by a 
combination of spatial inversion and an additional discrete symmetry. 



The same reasoning can be repeated for the case (eiy,?7vi/) = One 
finds that the space of projectors is the quaternionic Grassmannian GAr_,Ar++Ar_ (EI) = 
Sp(A^+ + A^_)/[Sp(A^+) X Sp(A^_)], whose homotopy group is again given in table 0. 

When {ewiVw) = ("l^+l); the Hamiltonian can be taken, in a suitable basis, 
real and skew symmetric, i.e., an element of so(n/). Any element of so(n/) can be 
transformed into a canonical form by an orthogonal transformation. For the Q-matrix, 
its canonical form Q is a matrix whose matrix elements are ±1, 

SQS^, 



Q 



where Q 



S e SO(n/), 
/ +1 
-1 







V 



(C.22) 



where we have assumed n/ is an even integer. We thus identified the space of projectors 
as SO(?7,/)/U(n//2). Finally, when {ew,Vw) = we identify the space of 

projectors as Sp(?7,j)/U(nj). 

To summarize, we have listed in table |C^, for each of the four cases of W- 
type symmetries, the space of projectors ^^G/H" together with their homotopy groups 
TidiG/H), which yields the classification of topological insulators (superconductors) 
protected by a combination of inversion and an additional discrete symmetry. Notice, we 
focus here on the implication of the combination of inversion and a discrete symmetries, 
but we do not require each symmetry separately. 



C.2.1. Example: {ewiflw) = (1? 1) in d = 2. We now specialize to the case 
of {ew,^w) = (l^l)- The relevant space of the projection operators is the real 
Grassmannian Gn_^_^n_^_+n_{M.) and the projection operator defines a map from the BZ 
onto the real Grassmannian. In particular, when d = 2, the projector Q{kx, ky) defines a 
map from S"^ or onto Gj\f^^N_^]\f_^_(R). The homotopy group Tr2[GN+,N-+N+(M)] = ^2 
tells us that the space of quantum ground states is partitioned into two topologically 
distinct sectors. 
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A representative Q-field configurations which is a non-trivial element of 
T^2[G N+,N-+N+{^)] = ^2 can be constructed following Ref. [^. For simplicity, take 
the case where we have four bands in total and two filled bands; take = = 2. 
Then such a representative Q-matrix is given by 





(0 





(0 
-ny 


(0 

-ny 
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ny 


(0 
-nx 






(0 




(0 

-ny 
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-ny 


(0 
-nV 





(0 


/ 



with 



n 



(0 



(^) := 



(0 „(0 ^(0 

X , ll'y , 



nl 



sin(/0) sin(6'), cos6' 



(C.23) 



cos(/0) sin(6'), 

where / G Z and 6 and </> is spherical coordinates parameterizing the BZ ~ S^. The 
projector is an non-trivial element of the second homotopy group when / is odd, while 
it is trivial when / is even. Observe that the normalized vector n^^\k) itself defines a 
map from S"^ onto S^, and wraps S"^ integer (= /) times. 

We can construct a Hamiltonian which has Q''^\k) as a projector. We consider the 
following four band tight-binding Hamiltonian on the square lattice: 

r 

-|-?/;^(r) hxip{r + x) + h.c. + t(j\r) hy t(j{r + y) + h.c] 
1^1, 0) and y = (0, 1), respectively, and 

Co \ / tao A(7y 



(C.24) 



where x 
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and t. A, /i G M are a parameter of the model. In /c-space. 
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where 



R{k) 



-2A sin ky, — 2A sin kx 



2t (cos kx + cos k^ 



(C.27) 



We will set t = A = 1 henceforth. As we change /i, there are four phases separated 
by quantum phase transitions at /i = ±4 and /i = 0. When > 4, the normalized 
vector R{k)/\R{k) \ does not wrap S"^ as we sweep the momentum and the system is in a 
trivial phase. On the other hand, for < 4 the normalized vector R{k)/\R{k)\ wraps 
the sphere 5*^ n = +1 times (or n = —1 times), and hence the system is in a non-trivial 
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phase. The projector constructed from ( |C.26| ) is topologically equivalent to Q^^\k) with 
/ odd. 

For topological insulators and superconductors protected by T, C, or combination 
of both, a diagnostics of bulk topological character is the existence of gapless edge 
modes. For insulators and superconductors with an inversion symmetry, it is less clear 
if looking for edge modes is useful to characterize the bulk topological character since 
boundary breaks inversion symmetry. [Nevertheless, as discussed recently in Ref. ||84|| , 
for the entanglement entropy spectrum, an inversion symmetry can protect the gapless 
entanglement entropy spectrum.] For the present case, however, the Hamiltonian has 
several accidental symmetries (see below), in particular T (odd), which may protect the 
gapless nature of edge states, if they exist, and if the number of Kramers pairs is odd. 
In Fig. pi] , the energy spectrum for = ±1 in the slab geometry is presented; here 
two edges along y-direction (located at x = and x = N^) are created, and the energy 
eigenvalues are plotted as a function of ky. Each eigenvalue is doubly degenerate for 
each ky. For fi = ±1, there are four edge modes, two of which localized at x = and 
the other two at x = N^,., while for > 4 there is no edge mode. 

The only symmetry we have assumed so far for ensembles of Hamiltonians is a 
W-tjpe symmetry with {ew,Vw) = (l^l)- The representative Hamiltonian ( |C.25 



C.26| ), however, accidentally has several other discrete symmetries. The Hamiltonian is 
invariant under the following discrete symmetries: Time-reversal symmetries T 

A~^n{-kyA = U{k), A = T,ao or T,ay, 
inversion symmetries X 

B~^ni-k)B = n{k), B = T,ao or r^a^, 
and charge- conjugation symmetries C 

Cni-kyC = -nik), C = T,a, or r,a,. 
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Although accidental, the presence of these symmetries allows us to interpret the 
Hamiltonian in many different ways. For example, since the Hamiltonian is invariant 
under T, A = TzfJy with = —A, it is a member of symmetry class AIL Furthermore, 
it supports two branches of modes per edge which are counter propagating. The 
representative Hamiltonian so constructed is thus a Z2 topological insulator in class 
All (i.e., it is in the QSH phase). With further imposing a C, C = T^cr^ with = +C, 
the single particle Hamiltonian can be interpreted as a member of symmetry class DHL 
Again, it is a non-trivial Z2 topological superconductor in class DHL Thus, non-trivial 
Z2 topological insulators (AH) and superconductors (Dili) in d = 2 dimensions can be 
interpreted, somewhat accidentally, a non-trivial element of vr2[G'Af_|_,Ar_+Af+(lR)] = Z2. 
To implement a W-tjpe symmetry, T for symmetry classes AH and Dili [A = r^ay 
with = -A) can be combined with an inversion symmetry represented by B = TzCFy. 
This inversion symmetry is, however, an artificial symmetry in the sense that it flips 
spin. 
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